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Table 6: Ages and metallicities for the best fit isochrones to the sample of 58 
double-lined eclipsing binaries for the different sets of tracks used. Systems 
YY Gem [la] and CM Dra [lb] are discussed in §4.2. T refers to the tracks 
used to derive the fit : (1) Geneva with overshooting, (2) Geneva without 
overshooting, (3) Padova, and (4) CG92. Maximum and minimum values 
are given at 68% c.l. (±1 a) for the fits of individual components, while the 
values for the combined system are at ±na. Ages consistent with the ZAMS 
are indicated by ’ZA’, that is, t = tzAMS in E<1- 1- 


# 

Name 


Component A 

Component B 

System A-l-B 

T 

r 


"^max 

^ ^min ■^max 

T 

^min 

^max 

z 

■^min ^max 

r 

"^min 

"^max 

Z 

■^min ■^max 

n 

2 

HS Aur 

1 

10.06 

9.95 

10.15 

0.040 0.020 0.040 

9.86 

9.24 

10.10 

0.040 

0.017 0.040 

9.98 

9.94 

10.02 

0.040 

0.035 0.040 

3 



2 

10.13 

9.93 

10.19 

0.040 0.020 0.040 

9.96 

9.28 

10.06 

0.039 

0.017 0.040 

9.98 

9.92 

10.02 

0.040 

0.035 0.040 

3 



3 

10.05 

7.43 

10.14 

0.049 0.022 0.100 

8.64 

ZA 

10.05 

0.073 

0.016 0.100 

9.98 

9.87 

10.06 

0.049 

0.021 0.050 

3 

3 

EW Ori 

1 

8.39 

ZA 

9.15 

0.040 0.024 0.040 

8.49 

ZA 

9.22 

0.040 

0.028 0.040 

8.42 

ZA 

8.79 

0.040 

0.036 0.040 




2 

7.35 

ZA 

9.11 

0.040 0.024 0.040 

7.27 

ZA 

9.15 

0.040 

0.025 0.040 

7.30 

ZA 

8.75 

0.040 

0.032 0.040 




3 

ZA 

ZA 

9.00 

0.051 0.023 0.069 

ZA 

ZA 

9.14 

0.053 

0.024 0.079 

ZA 

ZA 

8.80 

0.051 

0.032 0.065 




4 

5.91 

ZA 

9.10 

0.030 0.017 0.030 

5.91 

ZA 

9.16 

0.030 

0.018 0.030 

9.60 

ZA 

9.60 

0.030 

0.021 0.030 


4 

FL Lyr 

1 

9.36 

8.77 

9.55 

0.040 0.020 0.040 

9.43 

ZA 

9.75 

0.040 

0.021 0.040 

9.38 

9.19 

9.50 

0.040 

0.027 0.040 




2 

9.37 

8.67 

9.55 

0.040 0.021 0.040 

9.55 

ZA 

9.78 

0.040 

0.018 0.040 

9.40 

9.09 

9.50 

0.040 

0.028 0.040 




3 

7.73 

ZA 

9.50 

0.062 0.020 0.086 

3.79 

ZA 

9.93 

0.081 

0.029 0.100 

7.69 

ZA 

9.41 

0.065 

0.021 0.084 


5 

AI Phe 

1 

9.66 

9.64 

9.76 

0.009 0.009 0.028 

9.60 

9.55 

9.68 

0.008 

0.004 0.012 

9.66 

9.64 

9.67 

0.010 

0.008 0.011 




2 

9.71 

9.60 

9.73 

0.025 0.008 0.040 

9.60 

9.59 

9.66 

0.006 

0.006 0.013 

9.64 

9.61 

9.65 

0.009 

0.008 0.013 




3 

9.68 

9.67 

9.70 

0.027 0.022 0.035 

9.60 

9.54 

9.70 

0.007 0.006 0.013 

9.60 

9.59 

9.63 

0.008 

0.008 0.011 




4 

9.67 

9.64 

9.68 

0.011 0.010 0.015 

9.65 

9.63 

9.68 

0.010 

0.010 0.013 

9.66 

9.64 

9.67 

0.011 

0.010 0.012 


6 

UX Men 

1 

9.42 

8.90 

9.49 

0.033 0.021 0.060 

9.44 

8.80 

9.49 

0.030 

0.018 0.062 

9.45 

8.95 

9.56 

0.031 

0.022 0.056 




2 

9.43 

9.15 

9.55 

0.031 0.021 0.040 

9.47 

9.10 

9.58 

0.027 0.019 0.040 

9.45 

9.29 

9.52 

0.032 

0.024 0.040 




3 

8.83 

6.86 

9.20 

0.060 0.024 0.092 

7.97 

7.33 

9.33 

0.069 

0.020 0.087 

9.26 

9.15 

9.38 

0.034 0.025 0.044 




4 

9.38 

9.08 

9.57 

0.030 0.017 0.030 

9.44 

9.12 

9.60 

0.024 0.016 0.030 

9.40 

9.27 

9.51 

0.027 0.020 0.030 


7 

V1143 Cyg 

1 

ZA 

ZA 

8.97 

0.040 0.028 0.042 

8.49 

ZA 

9.12 

0.040 

0.027 0.045 

7.32 

ZA 

8.50 

0.040 

0.036 0.042 

1 



3 

ZA 

ZA 

8.84 

0.051 0.032 0.065 

ZA 

ZA 

8.84 

0.048 

0.029 0.064 

6.00 

ZA 

8.28 

0.050 

0.043 0.054 

1 



4 

5.40 

ZA 

8.88 

0.030 0.018 0.030 

7.92 

ZA 

9.00 

0.030 

0.018 0.030 

8.90 

ZA 

8.90 

0.030 

0.020 0.030 

2 

8 

DM Vir 

1 

9.26 

9.07 

9.34 

0.031 0.022 0.060 

9.26 

9.11 

9.34 

0.031 

0.020 0.060 

9.24 

9.05 

9.33 

0.034 0.023 0.055 

P 



3 

9.00 

8.80 

9.27 

0.055 0.026 0.062 

9.00 

8.80 

9.28 

0.055 

0.024 0.064 

9.15 

8.85 

9.25 

0.041 

0.028 0.063 

P 



4 

9.27 

9.04 

9.39 

0.027 0.017 0.030 

9.27 

9.02 

9.40 

0.027 0.016 0.030 

9.24 

9.15 

9.37 

0.030 

0.019 0.030 

n 

9 

BW Aqr 

1 

9.33 

9.29 

9.39 

0.026 0.019 0.037 

9.37 

9.29 

9.44 

0.017 0.013 0.025 

9.37 

9.31 

9.42 

0.021 

0.017 0.030 

1 



3 

9.10 

9.04 

9.40 

0.054 0.010 0.067 

9.34 

6.63 

9.38 

0.018 

0.003 0.100 

9.31 

9.24 

9.36 

0.024 0.018 0.036 

2 



4 

9.36 

9.24 

9.41 

0.023 0.017 0.030 

9.40 

9.32 

9.49 

0.016 

0.012 0.030 

9.38 

9.34 

9.42 

0.019 

0.016 0.025 

2 

10 

RZ Cha 

1 

9.32 

9.26 

9.35 

0.016 0.011 0.023 

9.33 

9.27 

9.35 

0.015 

0.010 0.022 

9.32 

9.29 

9.34 

0.015 

0.012 0.019 

P 



3 

9.30 

8.96 

9.31 

0.019 0.014 0.083 

9.30 

8.98 

9.33 

0.018 

0.015 0.088 

9.30 

9.27 

9.34 

0.018 

0.015 0.025 

P 



4 

9.34 

9.27 

9.35 

0.01^ 0.011 0.030 

9.35 

9.27 

9.38 

0.017 0.011 0.026 

9.34 

9.31 

9.35 

0.017 0.013 0.019 

n 

11 

V442 Cyg 

1 

9.21 

9.12 

9.26 

0.016 0.010 0.024 

9.24 

9.09 

9.36 

0.015 

0.011 0.023 

9.21 

9.17 

9.23 

0.015 

0.014 0.017 

p 



3 

8.88 

8.79 

9.23 

0.071 0.014 0.085 

8.60 

6.86 

9.28 

0.080 

0.012 0.100 

9.17 

9.14 

9.21 

0.018 

0.015 0.020 

p 



4 

9.23 

9.17 

9.26 

0.015 0.012 0.019 

9.26 

9.11 

9.34 

0.014 0.011 0.018 

9.23 

9.21 

9.26 

0.015 

0.013 0.016 

n 

12 

PV Pup 

1 

6.31 

ZA 

9.08 

0.040 0.019 0.058 

ZA 

ZA 

9.08 

0.040 

0.018 0.058 

7.53 

ZA 

8.96 

0.040 

0.024 0.053 

p 



3 

ZA 

ZA 

9.03 

0.050 0.022 0.091 

ZA 

ZA 

9.03 

0.050 

0.020 0.087 

5.00 

5.00 

8.80 

0.040 

0.028 0.072 

p 



4 

5.80 

ZA 

9.12 

0.030 0.017 0.030 

4.20 

ZA 

9.08 

0.030 

0.016 0.030 

8.50 

ZA 

8.75 

0.030 

0.021 0.030 

n 

13 

MY Cyg 

1 

8.95 

8.85 

9.05 

0.040 0.029 0.042 

8.95 

8.85 

9.04 

0.040 

0.034 0.043 

8.95 

8.85 

9.00 

0.040 

0.037 0.043 

1 



3 

8.91 

ZA 

9.05 

0.049 0.030 0.100 

8.89 

8.66 

9.08 

0.049 

0.037 0.054 

8.89 

8.83 

8.94 

0.050 

0.047 0.051 

2 



4 

9.05 

8.89 

9.16 

0.030 0.017 0.030 

9.03 

8.86 

9.16 

0.030 

0.018 0.030 

9.03 

9.02 

9.04 

0.029 

0.029 0.030 

2 

14 

RS Cha 

1 

8.92 

8.75 

8.99 

0.015 0.009 0.024 

9.01 

8.94 

9.02 

0.013 

0.008 0.023 

8.97 

8.93 

9.00 

0.015 

0.013 0.016 

1 



3 

8.89 

6.00 

9.00 

0.016 0.010 0.100 

8.69 

8.55 

9.00 

0.072 

0.010 0.089 

8.93 

8.82 

8.98 

0.016 

0.011 0.024 

2 



4 

8.95 

8.83 

9.01 

0.015 0.010 0.022 

9.04 

8.91 

9.08 

0.014 0.010 0.020 

8.98 

8.93 

9.03 

0.015 

0.011 0.019 

2 
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Table ^ : (continued) 


# 

Name 

T 

Component A 

Component B 

System A+B 

r 

^min 

^max 

Z 

■^min ^max 

T 

^min 

"^max 

Z Zmin 

■^max 

T 


^max 

^ ^min ^max 

n 

16 

WW Aur 

1 

8.31 

ZA 

8.69 

0.040 

0.022 0.060 

8.30 

ZA 

8.99 

0.040 0.014 0.060 

8.46 

ZA 

8.59 

0.039 0.030 0.058 

1 



3 

7.71 

ZA 

8.73 

0.050 

0.029 0.092 

8.28 

ZA 

8.97 

0.056 0.018 

0.100 

8.03 

ZA 

8.64 

0.050 0.032 0.078 

1 



4 

8.61 

ZA 

8.93 

0.030 

0.018 0.030 

8.76 

7.18 

9.06 

0.030 0.015 

0.030 

8.65 

7.60 

8.90 

0.030 0.021 0.030 

2 

17 

EK Cep 

1 

ZA 

ZA 

8.34 

0.023 

0.017 0.038 













3 

ZA 

ZA 

8.19 

0.026 

0.018 0.042 













4 

ZA 

ZA 

8.44 

0.020 

0.016 0.030 











18 

TZ For 

1 

9.00 

8.95 

9.05 

0.009 

0.007 0.040 

9.06 

9.00 

9.06 

0.028 0.012 

0.040 

9.04 

9.01 

9.08 

0.015 0.014 0.019 

1 



3 

8.98 

8.91 

9.04 

0.041 

0.007 0.044 

8.73 

8.70 

9.10 

0.064 0.020 

0.070 

9.03 

8.99 

9.03 

0.038 0.021 0.042 

1 



4 






9.04 

9.02 

9.12 

0.016 0.014 0.030 






19 

VV Pyx 

1 

8.75 

8.69 

8.79 

0.008 

0.006 0.011 

8.75 

8.69 

8.79 

0.008 0.006 

0.011 

8.75 

8.71 

8.79 

0.008 0.006 0.010 

1 



3 

8.73 

6.74 

8.84 

0.008 

0.006 0.100 

8.73 

6.74 

8.84 

0.008 0.006 

0.100 

8.72 

8.66 

8.76 

0.008 0.006 0.010 

1 



4 

8.73 

8.66 

8.80 

0.010 

0.010 0.013 

8.73 

8.67 

8.80 

0.010 0.010 

0.013 

8.76 

8.72 

8.80 

0.010 0.010 0.011 

1 

20 

AI Hya 

1 

8.93 

8.91 

8.97 

0.040 

0.015 0.040 

8.94 

8.90 

9.00 

0.040 0.026 

0.040 

8.94 

8.91 

8.97 

0.040 0.033 0.040 

1 



3 

8.83 

8.80 

8.95 

0.050 

0.022 0.050 

8.78 

8.70 

8.95 

0.050 0.048 

0.050 

8.92 

8.92 

8.97 

0.050 0.033 0.050 

1 



4 

8.98 

8.93 

9.00 

0.020 

0.019 0.023 

9.01 

8.96 

9.07 

0.030 0.018 

0.030 

9.02 

8.99 

9.03 

0.030 0.027 0.030 

1 

21 

EE Peg 

1 

8.35 

ZA 

8.49 

0.040 

0.022 0.058 

8.55 

ZA 

9.44 

0.040 0.013 

0.070 

8.34 

ZA 

8.60 

0.040 0.027 0.052 

1 



3 

8.19 

ZA 

8.55 

0.051 

0.030 0.080 

ZA 

ZA 

9.20 

0.035 0.014 0.093 

8.28 

ZA 

8.60 

0.045 0.027 0.072 

1 



4 

8.50 

7.96 

8.80 

0.030 

0.018 0.030 

8.53 

ZA 

9.31 

0.030 0.012 

0.030 

8.52 

7.78 

8.75 

0.030 0.017 0.030 

2 

22 

V1647 Sgr 

1 

8.23 

ZA 

8.62 

0.019 

0.011 0.075 

8.23 

ZA 

8.68 

0.016 0.010 

0.080 

8.24 

ZA 

8.49 

0.017 0.014 0.023 

1 



3 

8.11 

ZA 

8.48 

0.020 

0.012 0.081 

8.02 

ZA 

8.77 

0.017 0.010 

0.082 

8.38 

ZA 

8.40 

0.014 0.012 0.033 

1 



4 

8.36 

ZA 

8.70 

0.017 0.013 0.030 

8.36 

ZA 

8.77 

0.015 0.010 

0.030 

8.36 

7.50 

8.62 

0.016 0.014 0.020 

1 

23 

GZ CMa 

1 

8.75 

8.49 

8.77 

0.016 

0.008 0.036 

8.75 

8.12 

8.83 

0.018 0.011 

0.040 

8.73 

8.61 

8.76 

0.018 0.011 0.026 

1 



3 

8.73 

ZA 

8.77 

0.017 0.011 0.100 

6.33 

ZA 

8.77 

0.089 0.012 

0.100 

8.72 

8.53 

8.78 

0.018 0.013 0.035 

1 



4 

8.80 

8.50 

8.65 

0.016 

0.010 0.030 

8.81 

8.23 

8.93 

0.017 0.011 

0.030 

8.80 

8.63 

8.85 

0.017 0.013 0.030 

1 

24 

V624 Her 

1 

8.74 

8.65 

8.81 

0.033 

0.020 0.040 

8.91 

8.76 

8.98 

0.017 0.012 

0.028 

8.77 

8.73 

8.81 

0.025 0.021 0.032 

2 



3 

8.67 

8.51 

8.78 

0.050 

0.021 0.053 

8.58 

5.53 

8.90 

0.065 0.017 0.100 

8.74 

8.68 

8.78 

0.030 0.021 0.042 

2 



4 

8.79 

8.75 

8.88 

0.030 

0.020 0.030 

8.96 

8.80 

9.00 

0.017 0.013 

0.030 

8.80 

8.78 

8.83 

0.028 0.024 0.030 

2 

25 

YZ Cas 

1 

8.66 

7.80 

8.68 

0.017 0.012 0.100 

8.88 

ZA 

9.37 

0.020 0.010 

0.100 

8.65 

8.48 

8.75 

0.018 0.014 0.027 

1 



3 

8.35 

7.94 

8.74 

0.062 

0.011 0.079 

8.60 

ZA 

9.27 

0.020 0.011 

0.100 

8.62 

8.43 

8.71 

0.020 0.013 0.036 

1 



4 

8.72 

8.42 

8.77 

0.017 0.011 0.030 

8.89 

ZA 

9.20 

0.017 0.010 

0.030 

8.73 

8.54 

8.80 

0.017 0.013 0.030 

1 

26 

SZ Cen 

1 

8.81 

8.58 

8.88 

0.013 

0.006 0.050 

8.83 

8.50 

8.90 

0.019 0.011 

0.050 

8.82 

8.80 

8.86 

0.049 0.017 0.052 

1 



3 

8.56 

8.48 

8.85 

0.056 

0.017 0.060 

8.56 

8.51 

8.85 

0.054 0.010 

0.058 

8.56 

8.51 

8.57 

0.055 0.053 0.060 

1 



4 

8.83 

8.81 

8.87 

0.014 0.011 0.030 

8.87 

8.77 

8.90 

0.024 0.016 

0.030 

8.87 

8.84 

8.87 

0.027 0.021 0.030 

1 

27 

f3 Aur 

1 

8.68 

8.59 

8.74 

0.013 

0.010 0.021 

8.70 

8.62 

8.74 

0.014 0.009 

0.021 

8.69 

8.63 

8.73 

0.014 0.010 0.020 

1 



3 

8.37 

7.63 

8.73 

0.066 

0.005 0.095 

8.38 

7.42 

8.75 

0.066 0.006 

0.094 

8.36 

8.05 

8.74 

0.066 0.007 0.076 

1 



4 

8.72 

8.65 

8.75 

0.015 

0.010 0.022 

8.74 

8.67 

8.77 

0.015 0.010 

0.024 

8.71 

8.68 

8.75 

0.017 0.012 0.019 

1 

28 

V1031 Ori 

1 

8.74 

8.40 

8.75 

0.019 ft013 0.070 

8.76 

8.25 

8.85 

0.022 0.008 

0.090 

8.77 

8.40 

8.80 

0.021 0.016 0.065 

1 



3 

8.65 

8.40 

8.75 

0.042 

0.007 0.060 

8.75 

8.30 

8.83 

0.021 0.009 

0.076 

8.75 

8.40 

8.78 

0.020 0.018 0.060 

1 



4 

8.77 

8.69 

8.79 

0.028 

0.016 0.030 

8.79 

8.56 

8.89 

0.025 0.010 

0.030 

8.77 

8.70 

8.78 

0.027 0.020 0.030 

1 

29 

TZ Men 

1 

8.11 

ZA 

8.56 

0.020 

0.011 0.082 

ZA 

ZA 

9.09 

0.019 0.011 

0.078 

8.15 

ZA 

8.49 

0.019 0.016 0.034 

1 



3 

8.03 

ZA 

8.41 

0.022 

0.011 0.082 

ZA 

ZA 

9.06 

0.021 0.013 

0.083 

8.17 

ZA 

8.46 

0.020 0.013 0.073 

1 



4 

8.24 

ZA 

8.55 

0.019 

0.012 0.030 

8.29 

ZA 

9.13 

0.017 0.011 

0.030 

8.25 

ZA 

8.50 

0.018 0.014 0.030 

1 

30 

WX Cep 

1 

8.69 

8.60 

8.75 

0.028 

0.016 0.040 

8.71 

8.25 

8.80 

0.019 0.010 

0.075 

8.72 

8.67 

8.74 

0.020 0.018 0.027 

1 



3 

8.65 

8.36 

8.76 

0.042 

0.016 0.057 

8.59 

8.31 

8.74 

0.035 0.011 

0.076 

8.69 

8.66 

8.72 

0.021 0.017 0.032 

1 



4 

8.73 

8.69 

8.78 

0.020 

0.017 0.030 

8.77 

8.56 

8.80 

0.019 0.011 

0.030 

8.76 

8.70 

8.76 

0.020 0.017 0.030 

1 

31 

V451 Oph 

1 

8.40 

ZA 

8.49 

0.016 

0.005 0.100 

8.20 

ZA 

8.68 

0.024 0.010 

0.090 

8.37 

ZA 

8.53 

0.020 0.012 0.100 

1 



3 

8.01 

ZA 

8.57 

0.063 

0.005 0.100 

8.07 

ZA 

8.65 

0.029 0.010 

0.100 

8.28 

ZA 

8.52 

0.024 0.011 0.092 

1 



4 

8.45 

ZA 

8.45 

0.016 

0.010 0.030 

8.30 

ZA 

8.70 

0.024 0.012 

0.030 

8.38 

7.10 

8.58 

0.021 0.014 0.030 

1 





























Table ^ : (continued) 


# 

Name 

T 

Component A 

Component B 

System A+B 

r 

^min 

^max 

Z 

^min ^max 

T 

^min 

"^max 

Z .^min 

■^max 

r 

^min 

^max 

Z 

^min ^max 

n 

32 

PV Cas 

1 

7.66 

ZA 

8.39 

0.040 

0.027 0.041 

7.94 

ZA 

8.41 

0.040 0.025 

0.042 

7.88 

ZA 

8.15 

0.040 

0.038 0.041 

1 



3 

7.74 

ZA 

8.19 

0.050 

0.032 0.055 

7.90 

ZA 

8.48 

0.050 0.029 

0.057 

3.68 

ZA 

7.87 

0.050 

0.046 0.051 

1 



4 

7.22 

ZA 

8.46 

0.030 

0.019 0.030 

7.89 

ZA 

8.53 

0.030 0.018 

0.030 

7.64 

ZA 

8.12 

0.030 

0.027 0.030 

1 

33 

IQ Per 

1 

ZA 

ZA 

7.88 

0.040 

0.025 0.040 

ZA 

ZA 

8.60 

0.032 0.021 

0.043 

ZA 

ZA 

8.00 

0.039 

0.020 0.041 

2 



3 

ZA 

ZA 

7.81 

0.049 

0.025 0.055 

ZA 

ZA 

8.44 

0.036 0.023 

0.060 

7.95 

ZA 

8.10 

0.020 

0.018 0.060 

2 



4 

6.23 

ZA 

8.00 

0.030 

0.018 0.030 

ZA 

ZA 

8.63 

0.030 0.018 

0.030 

6.20 

ZA 

7.80 

0.030 

0.024 0.030 

3 

34 

Hya 

1 

8.27 

7.99 

8.30 

0.013 

0.008 0.021 

8.26 

7.92 

8.49 

0.015 0.011 

0.020 

8.27 

8.24 

8.30 

0.014 0.011 0.018 

1 



3 

7.92 

7.90 

8.30 

0.059 

0.008 0.063 

8.24 

ZA 

8.32 

0.015 0.010 

0.087 

8.28 

8.27 

8.30 

0.014 0.012 0.017 

1 



4 

8.31 

8.21 

8.33 

0.015 

0.010 0.030 

8.32 

7.85 

8.50 

0.016 0.011 

0.030 

8.31 

8.29 

8.33 

0.015 

0.013 0.017 

1 

35 

C Phe 

1 

7.92 

7.73 

8.04 

0.011 

0.007 0.019 

8.21 

7.85 

8.39 

0.007 0.005 

0.008 

7.98 

7.88 

8.04 

0.008 

0.007 0.009 

2 



3 

6.89 

ZA 

8.08 

0.070 

0.007 0.084 

8.17 

7.75 

8.33 

0.007 0.005 

0.009 

8.01 

7.98 

8.05 

0.008 

0.007 0.008 

1 



4 

7.99 

7.80 

8.09 

0.013 

0.010 0.019 

8.17 

7.69 

8.35 

0.010 0.010 

0.011 

8.03 

7.85 

8.11 

0.010 

0.010 0.012 

3 

36 

GG Lup 

1 

ZA 

ZA 

7.64 

0.025 

0.014 0.061 

ZA 

ZA 

8.43 

0.016 0.009 

0.100 

ZA 

ZA 

7.50 

0.023 

0.015 0.037 

1 



3 

ZA 

ZA 

7.66 

0.021 

0.013 0.049 

ZA 

ZA 

8.40 

0.018 0.008 

0.077 

ZA 

ZA 

7.58 

0.020 

0.014 0.035 

1 



4 

ZA 

ZA 

7.67 

0.022 

0.014 0.030 

ZA 

ZA 

8.40 

0.016 0.010 

0.030 

6.45 

ZA 

7.70 

0.020 

0.015 0.030 

1 

37 

V760 Sco 

1 

7.51 

ZA 

7.77 

0.015 

0.006 0.068 

6.93 

ZA 

7.74 

0.020 0.008 

0.056 

7.36 

6.60 

7.59 

0.017 0.012 0.027 

1 



3 

7.59 

ZA 

7.75 

0.011 

0.006 0.074 

7.20 

ZA 

7.75 

0.015 0.011 

0.065 

7.52 

7.10 

7.68 

0.013 

0.009 0.019 

1 



4 

7.62 

ZA 

7.68 

0.016 

0.010 0.030 

7.34 

ZA 

7.68 

0.018 0.011 

0.030 

7.56 

7.10 

7.71 

0.017 0.013 0.030 

1 

38 

DI Her 

1 

ZA 

ZA 

7.60 

0.025 

0.009 0.050 

ZA 

ZA 

7.62 

0.040 0.014 0.048 

ZA 

ZA 

7.41 

0.033 

0.016 0.048 

1 



3 

ZA 

ZA 

7.56 

0.021 

0.009 0.060 

ZA 

ZA 

7.67 

0.028 0.013 

0.058 

5.00 

ZA 

7.49 

0.023 

0.013 0.054 

1 



4 

ZA 

ZA 

7.64 

0.022 

0.011 0.030 

ZA 

ZA 

7.64 

0.020 0.014 0.030 

ZA 

ZA 

7.56 

0.020 

0.014 0.030 

1 

39 

U Oph 

1 

7.55 

ZA 

7.87 

0.023 

0.007 0.070 

7.47 

ZA 

7.87 

0.023 0.007 0.070 

7.57 

ZA 

7.70 

0.021 

0.012 0.072 

1 



3 

7.17 

ZA 

7.88 

0.053 

0.005 0.073 

7.60 

ZA 

7.93 

0.017 0.006 

0.073 

7.66 

ZA 

7.87 

0.016 

0.007 0.070 

1 



4 

7.56 

ZA 

7.61 

0.030 

0.011 0.030 

7.49 

ZA 

7.97 

0.030 0.010 

0.030 

7.59 

6.88 

7.81 

0.026 

0.010 0.030 

1 

40 

GV Vel 

1 

7.42 

6.80 

7.56 

0.031 

0.012 0.054 

7.39 

6.80 

7.55 

0.033 0.012 

0.054 

7.42 

7.10 

7.59 

0.031 

0.015 0.050 

1 



3 

7.57 

6.86 

7.71 

0.019 

0.010 0.058 

7.56 

6.64 

7.71 

0.020 0.010 

0.058 

7.52 

7.00 

7.65 

0.023 

0.012 0.053 

1 



4 

7.62 

7.33 

7.71 

0.021 

0.013 0.030 

7.62 

7.29 

7.71 

0.021 0.014 0.030 

7.64 

7.45 

7.68 

0.020 

0.017 0.030 

1 

41 

V539 Ara 

1 

7.59 

7.00 

7.70 

0.017 0.007 0.062 

7.57 

ZA 

7.78 

0.017 0.007 0.068 

7.58 

7.10 

7.70 

0.017 0.010 0.060 

1 



3 

7.66 

7.05 

7.76 

0.012 

0.006 0.063 

6.95 

ZA 

7.82 

0.055 0.006 

0.071 

7.62 

7.05 

7.72 

0.015 

0.004 0.061 

1 



4 

7.59 

7.44 

7.74 

0.027 0.010 0.030 

7.66 

6.72 

7.82 

0.020 0.010 

0.030 

7.65 

7.44 

7.73 

0.020 

0.012 0.030 

1 

42 

QX Gar 

1 

6.17 

ZA 

7.18 

0.040 

0.008 0.055 

6.42 

ZA 

7.18 

0.036 0.010 

0.052 

6.16 

ZA 

7.10 

0.040 

0.014 0.070 

1 



3 

6.70 

ZA 

7.11 

0.022 

0.007 0.059 

6.49 

ZA 

7.22 

0.030 0.008 

0.059 

5.89 

ZA 

7.12 

0.040 

0.010 0.054 

1 



4 

6.92 

ZA 

7.19 

0.024 0.011 0.030 

7.03 

ZA 

7.32 

0.022 0.012 

0.030 

7.00 

6.20 

7.13 

0.020 

0.016 0.030 

1 

43 

GW Cep 

1 

6.71 

ZA 

6.99 

0.020 

0.003 0.065 

6.95 

ZA 

7.10 

0.008 0.002 

0.080 

6.81 

ZA 

7.00 

0.014 0.006 0.071 

1 



3 

6.73 

ZA 

7.03 

0.018 

0.004 0.063 

6.98 

ZA 

7.03 

0.007 0.002 

0.074 

6.26 

ZA 

7.00 

0.040 

0.006 0.070 

1 



4 

6.84 

5.80 

6.96 

0.020 

0.010 0.030 

6.96 

ZA 

6.96 

0.013 0.010 

0.030 

6.84 

6.36 

7.00 

0.022 

0.010 0.030 

1 

44 

V478 Cyg 

1 

6.91 

6.00 

7.00 

0.007 otiOl 0.070 

6.95 

6.00 

6.99 

0.004 0.001 

0.078 

6.93 

6.00 

7.00 

0.006 

0.001 0.080 

1 



3 

6.49 

5.87 

6.96 

0.046 

0.004 0.058 

6.45 

6.01 

6.96 

0.052 0.004 0.060 

6.46 

6.20 

7.00 

0.050 

0.004 0.054 

1 



4 

6.78 

6.55 

6.99 

0.027 0.010 0.030 

6.78 

6.55 

6.99 

0.027 0.010 

0.030 

6.76 

6.65 

6.90 

0.030 

0.010 0.030 

1 

45 

EM Car 

1 

6.80 

ZA 

6.89 

0.004 0.001 0.100 

6.79 

ZA 

6.90 

0.004 0.001 

0.100 

6.81 

5.00 

6.84 

0.003 

0.001 0.090 

1 



3 

6.29 

ZA 

6.89 

0.052 

0.000 0.062 

6.23 

ZA 

6.94 

0.052 0.000 

0.064 

6.79 

5.60 

6.88 

0.004 0.001 0.057 

1 



4 

6.61 

6.26 

6.84 

0.030 

0.010 0.030 

6.61 

5.92 

6.84 

0.027 0.010 

0.030 

6.60 

6.40 

6.82 

0.030 

0.010 0.030 

1 

46 

CG Cyg 

1 

9.46 

ZA 

9.81 

0.040 

0.012 0.058 

9.63 

8.21 

9.90 

0.040 0.019 

0.042 

9.53 

8.80 

9.68 

0.040 

0.025 0.041 

1 



2 

9.53 

ZA 

9.83 

0.040 

0.012 0.040 

9.51 

8.17 

9.89 

0.040 0.018 

0.040 

9.50 

8.71 

9.78 

0.040 

0.021 0.040 

1 



3 

8.24 

ZA 

9.81 

0.059 

0.012 0.097 

5.00 

ZA 

10.14 

0.066 0.034 0.084 

8.00 

ZA 

9.76 

0.064 0.036 0.082 

1 




























Table ^ : (continued) 


# 

Name 

T 

Component A 

Component B 

System A+B 

r 

^min 

^max 

Z 

^min ^max 

T 

^min 

"^max 

Z Zjnin 

■^max 

T 


^max 

^ ^min ^max 

n 

47 

RT And 

1 

9.37 

ZA 

9.66 

0.040 

0.013 0.060 

ZA 

ZA 

9.40 

0.040 0.022 

0.040 

ZA 

ZA 

9.20 

0.040 0.036 0.045 

1 



2 

9.30 

ZA 

9.60 

0.040 

0.014 0.040 

ZA 

ZA 

9.39 

0.040 0.021 

0.040 

7.18 

ZA 

9.40 

0.040 0.023 0.040 

1 



3 

7.71 

ZA 

9.49 

0.061 

0.018 0.100 

7.50 

ZA 

9.33 

0.051 0.024 0.065 

7.00 

ZA 

9.31 

0.051 0.026 0.066 

1 

48 

AR Aur 

1 

7.28 

ZA 

8.43 

0.016 

0.007 0.031 

8.21 

ZA 

8.65 

0.013 0.006 

0.031 

7.90 

ZA 

8.75 

0.014 0.011 0.018 

1 



3 

6.63 

ZA 

8.33 

0.016 

0.009 0.032 

8.13 

ZA 

8.41 

0.014 0.007 0.035 

7.78 

ZA 

8.20 

0.016 0.012 0.019 

1 



4 

7.77 

ZA 

8.41 

0.015 

0.010 0.030 

8.28 

ZA 

8.48 

0.013 0.010 

0.030 

8.07 

7.10 

8.30 

0.014 0.012 0.016 

1 

49 

AG Per 

1 

7.48 

ZA 

7.77 

0.010 

0.003 0.088 

7.12 

ZA 

7.77 

0.012 0.004 0.065 

7.38 

ZA 

7.67 

0.010 0.006 0.021 

1 



3 

7.55 

ZA 

7.57 

0.008 

0.000 0.090 

7.31 

ZA 

7.64 

0.009 0.000 

0.078 

7.49 

ZA 

7.71 

0.008 0.005 0.066 

1 



4 

7.54 

ZA 

7.61 

0.012 

0.010 0.030 

7.31 

ZA 

7.68 

0.013 0.010 

0.030 

7.45 

6.40 

7.70 

0.013 0.010 0.019 

1 

50 

YCyg 

1 

6.76 

ZA 

6.91 

0.002 

0.001 0.015 

6.75 

ZA 

6.91 

0.002 0.001 

0.014 

6.75 

5.04 

6.90 

0.002 0.001 0.007 

1 



3 

6.77 

ZA 

6.98 

0.002 

0.000 0.075 

6.76 

ZA 

6.98 

0.002 0.000 

0.077 

6.79 

ZA 

6.97 

0.002 0.000 0.071 

1 



4 

5.97 

ZA 

6.60 

0.010 

0.010 0.016 

6.05 

ZA 

6.60 

0.010 0.010 

0.016 

6.00 

ZA 

6.47 

0.010 0.010 0.013 

2 

51 

AH Cep 

1 

6.47 

ZA 

6.99 

0.037 0.002 0.070 

6.82 

ZA 

6.99 

0.015 0.002 

0.090 

6.49 

ZA 

7.00 

0.036 0.004 0.068 

1 



3 

6.55 

ZA 

6.88 

0.025 

0.004 0.058 

6.41 

ZA 

6.88 

0.038 0.004 0.064 

6.46 

ZA 

7.00 

0.033 0.004 0.055 

1 



4 

6.82 

5.29 

6.96 

0.020 

0.010 0.030 

6.85 

ZA 

6.93 

0.021 0.010 

0.030 

6.83 

6.34 

7.00 

0.020 0.010 0.030 

1 

52 

V3903 Sgr 

1 

6.50 

ZA 

6.72 

0.007 0.001 0.061 

6.50 

ZA 

6.80 

0.009 0.001 

0.052 

6.40 

ZA 

6.76 

0.010 0.001 0.068 

1 



3 

6.52 

ZA 

6.78 

0.004 0.000 0.064 

5.29 

ZA 

7.02 

0.062 0.000 

0.083 

5.00 

ZA 

6.75 

0.055 0.001 0.067 

1 



4 

6.07 

ZA 

6.61 

0.025 

0.010 0.030 

6.70 

5.71 

6.83 

0.010 0.010 

0.030 

6.52 

5.70 

6.68 

0.010 0.010 0.030 

1 

53 

DH Cep 

1 

6.49 

ZA 

6.65 

0.003 

0.001 0.100 

6.57 

ZA 

6.72 

0.003 0.001 

0.100 

6.51 

ZA 

6.66 

0.002 0.001 0.100 

1 



3 

6.48 

ZA 

6.52 

0.002 

0.000 0.064 

6.55 

ZA 

6.52 

0.003 0.000 

0.066 

ZA 

ZA 

6.65 

0.053 0.000 0.058 

1 



4 

6.42 

ZA 

6.60 

0.010 

0.010 0.030 

6.50 

ZA 

6.65 

0.010 0.010 

0.030 

6.46 

4.44 

6.56 

0.010 0.010 0.030 

1 

54 

V477 Cyg 

1 

8.49 

ZA 

9.05 

0.012 

0.004 0.036 

ZA 

ZA 

9.43 

0.021 0.007 0.058 

ZA 

ZA 

8.79 

0.021 0.008 0.035 

1 



3 

7.58 

ZA 

9.12 

0.017 0.003 0.100 

8.33 

ZA 

9.60 

0.017 0.001 

0.100 

8.12 

ZA 

9.18 

0.018 0.002 0.100 

1 



4 

8.57 

ZA 

9.00 

0.012 

0.010 0.030 

7.35 

ZA 

9.37 

0.017 0.010 

0.030 

7.98 

ZA 

8.93 

0.018 0.010 0.030 

1 

55 

YY Sgr 

1 

7.77 

ZA 

8.13 

0.011 

0.004 0.100 

7.88 

ZA 

8.12 

0.009 0.004 0.100 

7.75 

ZA 

8.10 

0.011 0.005 0.027 

1 



3 

8.01 

ZA 

8.20 

0.004 0.001 0.100 

8.10 

ZA 

8.28 

0.004 0.001 

0.100 

7.81 

ZA 

8.18 

0.010 0.003 0.100 

1 



4 

7.86 

ZA 

8.14 

0.013 

0.010 0.030 

7.95 

ZA 

8.23 

0.011 0.010 

0.030 

7.83 

ZA 

8.11 

0.013 0.010 0.026 

1 

56 

V526 Sgr 

1 

8.36 

ZA 

8.63 

0.013 

0.007 0.025 

8.82 

8.32 

8.98 

0.007 0.004 0.011 

8.52 

8.41 

8.62 

0.010 0.009 0.011 

1 



3 

8.29 

ZA 

8.77 

0.014 0.002 0.100 

8.74 

ZA 

9.08 

0.007 0.002 

0.100 

8.44 

ZA 

8.75 

0.011 0.005 0.100 

1 



4 

8.43 

7.50 

8.65 

0.014 0.010 0.020 

8.72 

ZA 

8.93 

0.010 0.010 

0.013 

8.59 

8.15 

8.77 

0.010 0.010 0.013 

1 

57 

AD Boo 

1 

9.24 

8.80 

9.36 

0.040 

0.031 0.042 

ZA 

ZA 

9.25 

0.040 0.025 

0.042 

9.18 

8.81 

9.34 

0.040 0.037 0.043 

2 



2 

9.21 

8.76 

9.35 

0.040 

0.029 0.040 

ZA 

ZA 

9.33 

0.040 0.026 

0.040 

9.16 

8.82 

9.31 

0.040 0.037 0.043 

2 



3 

ZA 

ZA 

6.47 

0.075 

0.059 0.093 

ZA 

ZA 

8.90 

0.057 0.041 

0.073 

ZA 

ZA 

9.65 

0.066 0.008 0.100 

1 

58 

SW CMa 

1 

8.80 

7.90 

8.90 

0.024 0.010 0.080 

8.83 

ZA 

9.00 

0.024 0.010 

0.087 

8.80 

8.20 

8.90 

0.024 0.015 0.085 

1 



3 

8.67 

ZA 

9.00 

0.086 

0.001 0.100 

8.81 

ZA 

9.00 

0.024 0.003 

0.100 

8.76 

ZA 

9.00 

0.031 0.003 0.100 

1 



4 

8.80 

8.60 

8.90 

0.030 

0.014 0.030 

8.85 

8.40 

9.02 

0.028 0.010 

0.030 

8.81 

8.66 

8.93 

0.030 0.015 0.030 

1 

59 

V505 Per 

1 

9.15 

ZA 

9.41 

0.019 (^012 0.088 

9.18 

ZA 

9.49 

0.023 0.012 

0.090 

9.11 

ZA 

9.40 

0.020 0.015 0.032 

1 



2 

9.14 

ZA 

9.39 

0.019 

0.013 0.033 

9.18 

ZA 

9.49 

0.023 0.013 

0.040 

9.11 

ZA 

9.38 

0.021 0.015 0.032 

1 



3 

8.93 

ZA 

9.61 

0.021 

0.003 0.100 

9.05 

ZA 

9.67 

0.021 0.003 

0.100 

8.95 

ZA 

9.60 

0.021 0.020 0.100 

1 



4 

9.12 

ZA 

9.42 

0.016 

0.012 0.026 

9.18 

ZA 

9.50 

0.018 0.012 

0.030 

9.10 

ZA 

9.40 

0.017 0.013 0.023 

1 
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Abstract. Detached, double-lined spectroscopic binaries which are also eclipsing provide the most accurate de¬ 
terminations of stellar mass, radius, temperature and distance-independent luminosity for each of their individual 
components, and hence constitute a stringent test of single-star stellar evolution theory. We compile a large sample 
of 60 non interacting, well-detached systems mostly with typical errors smaller than 2% for mass and radius and 
smaller than 5% for effective temperature, and compare them with the properties predicted by stellar evolution¬ 
ary tracks from a minimization method. To assess the systematic errors introduced by a given set of tracks, we 
compare the results obtained using three widely-used independent sets of tracks, computed with different physi¬ 
cal ingredients (the Geneva, Padova and Granada models). We also test the hypothesis that the components of 
these systems are coeval and have the same metallicity, and compare the derived ages and metallicities with the 
ones obtained by fitting a single isochrone to the system. Overall, there is a good agreement among the different 
determinations, and we provide a comprehensive discussion on the sub-sample of systems which either present 
problems or have estimated metallicities. Although within the errors the published tracks can fit most of the 
systems, a large degeneracy between age and metallicity remains. The power of the test is thus limited because 
the metallicities of most of the systems are unknown. 

Key words, stars: evolution - stars: fundamental parameters - stars: binaries - binaries: eclipsing - stars: HR 
diagram - Galaxy: evolution - stars: statistics 


1. Introduction 

The theory of stellar evolution is one of the most suc¬ 
cessful in astrophysics, and has been widely applied to 
phenomena ranging from stellar oscillations to the evolu¬ 
tion of galaxies. While open and globular clusters have 
been the classical tests of the theory because their colour- 
magnitude diagrams (CMD) constrain the shape of an 
isochrone in detail, they are not ideal probes, since, 
among other things, they are contaminated by held stars, 
have unresolved binaries, suffer dynamical evolution, may 
present differential extinction, and may be distant enough 
to prevent very accurate determinations of their dis¬ 
tances (see Andersen & Nordstrom 1999 for a recent re¬ 
view on the observational requirements for open clusters). 
Furthermore, Hipparcos measures show that there is a 
very fuzzy correlation, if any, between the main sequence 
position and metallicity (van Leeuwen 1999, Lebreton et 
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al. 1999, and references therein), casting some doubts - 
at least as far as the metallicity scaling is concerned- on 
current stellar evolution models. 

So far relatively less attention has been paid, in com¬ 
parison, to binary systems. Yet these systems provide the 
only way to measure directly stellar masses, and, if the 
components do not interact strongly and suffer mass trans¬ 
fer, they should have independent evolution and be rep¬ 
resentative of single stars. Visual systems are hampered 
by the lack of radial velocity measures (but see Pourbaix 
2000), and hence mass ratios are measured indirectly (e.g. 
Fernandes et al. 1998 for a few examples). Interferometric 
systems, although promising, are too few to provide a good 
sample (e.g. Hummel et al. 1995, 2001, Torres et al. 2002). 
In this context, double-lined eclipsing binaries provide the 
most accurate determinations of stellar radii, masses and 
temperatures, and hence are ideal stringent tests of the 
theory, as first pointed out by Stromgren (1967) (see also 
Andersen 1997). Testing evolutionary tracks is important 
not only for the understanding of the properties of the 
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stars, but also to assess the uncertainties in, for example, 
the interpretation of the spectro-photometric properties of 
galaxies, from low to high redshifts, or the chemical evo¬ 
lution of our Galaxy. 

Although only a handful of systems are extremely well- 
measured (see below), with typical errors in masses and 
radii of less than 2 per cent (Andersen 1991), the sam¬ 
ple of systems is bound to increase in the near future, 
with dedicated observing programs (at the Danish 50cm 
SAT, see Clausen, Kelt & Olsen 2001; Elodie at OHP, 
e.g. Kurpinska-Winiarska & Oblak 2000), the analysis of 
Hipparcos data (e.g., Fabricius et al. 2002), the systematic 
detections of thousands of systems as a by-product of mi- 
crolensing surveys (see, e.g., Alcock et al. 1997, and Ferlet, 
Maillard & Raban 1997 for an overview), and with future 
space missions such as DIVA, SIM and GAIA. In addition, 
these systems have accurate absolute dimensions and pro¬ 
vide another distance indicator, independent of any dis¬ 
tance calibration (e.g. Paczyhski 1997, Fitzpatrick et al. 
2002). However, in practice, the power of the test using 
these systems is limited by the small volume in parameter 
space sampled by most of the published theoretical tracks. 
For instance, many tracks adopt a single helium abun¬ 
dance Y for a given metallicity Z, and do not consider 
any possible dispersion, while others hx the overshooting 
parameter or the mixing length parameter to prescribed 
values. Even though these values provide good hts to the 
colour magnitude diagrams of clusters, it is likely that 
these parameters change with stellar mass or metal abun¬ 
dance. The extant (and growing) sample of binaries will 
yield extremely important constraints on these parame¬ 
ters (and their possible correlations) when tracks with a 
wider range of properties will be available. 

The aims of this paper are three-fold. First, we want to 
test whether current evolutionary tracks can account for 
the properties of these extremely well measured systems. 
This is a requisite before using them for further inferences. 
Even though a given set of tracks may seemingly give good 
results, and provide good fits, it will not produce an abso¬ 
lute determination because of the underlying hypotheses 
used in the computation of the tracks. For instance the 
determination of the Helium abundance relies critically 
on the effective temperature scale, and so different tracks, 
using different scales, may yield different abundances or 
ages. To assess the possible systematic effects introduced 
by a given set, we take three widely-used sets of tracks 
which have been computed independently and with dif¬ 
ferent physical ingredients. For the same data set, which 
includes only the best measured stars currently available, 
are these sets consistent with each other? 

Secondly, we want to test to which extent stellar evo¬ 
lution is a predictive theory, in the sense that given some 
observables it can predict the remaining ones, if the predic¬ 
tions are unique. This relies on the so-called Vogt-Russell 
theorem^ (Vogt 1926) which states that stars of the same 

^ We did not find any convincing reason to add Russell’s 
name to Vogt’s theorem, given that even though Russell may 


chemical composition (and hence same opacity and energy 
generation rate) have radii, effective temperatures and lu¬ 
minosities solely determined by their masses. Although 
this theorem has been widely criticized, mainly because 
counterexamples were found (albeit in rather contrived 
situations, see e.g. Lauterborn 1972, 1973; Kahler 1978), 
it remains the basis of stellar evolution theory and is im- 
plicitely assumed in all analyses of stellar evolution data. 
As far as we know, no tests have been carried out to assess 
empirically the validity of this theorem. So, given, say, a 
luminosity and effective temperature, can the tracks pro¬ 
vide a unique value for the mass and the radius? Are these 
values consistent with the measured ones? Allende Prieto 
& Lambert (1999) found that one set of tracks was able 
to reproduce, assuming solar metallicity, radii and masses 
with errors smaller than 8%, and effective temperatures 
within 2% for mainly main sequence stars. In a different 
study. Young et al. (2001) found agreement with the mea¬ 
sured radii within 3%, and temperatures within 4% but 
luminosities could only be reproduced to 11% with their 
set of tracks. Are these figures the same for other pub¬ 
lished tracks or for stars with lower gravities? 

Furthermore, and this is the third and main goal 
of the paper, can these sets of tracks predict the ages 
and metallicities of these systems? Proper testing of the 
accuracy of stellar dating has not yet been done before in 
detail (see e.g. von Hippel et al., 2001 for a summary). 
The components of these detached systems provide an 
ideal sample not only to properly calibrate the stellar 
tracks, but also to assess the uncertainties in the stellar 
ages used for the study of the evolution of our Galaxy. 
For instance Edvardsson et al. (1993) used a set of tracks 
from VandenBerg (1985) with a fixed set of physical 
parameters, and succeeded in reaching an accuracy of 
about 0.1 dex in the relative dating of F dwarfs. We want 
to assess whether this accuracy is robust to changes in the 
evolutionary tracks. In this context it is curious to observe 
that there has been no systematic test to see whether the 
components of these systems are indeed coeval. Popper 
(1997) observed that some secondary stars between 0.7 
and 1.1 Mq appeared systematically older (by at least a 
factor of two) than their primaries. Are the ages derived 
for each component consistent with them being on a single 
isochrone? Are their metallicities the same? What are the 
uncertainties (including systematics) in the absolute ages? 

This paper presents the sample, the method and the 
detailed hts, discusses the implications of the results on 
the coevality of the components, the ages and metallici¬ 
ties inferred, and their relation to the chemical evolution 
of the Galaxy. When relevant we also present comparisons 
with the results obtained by Pols et al. (1997), Ribas et 
al. (2000), and Young et al. (2001), hereafter referred to 
as P97, ROO and YOl respectively. Section §2 briefly de- 


had derived it independently of Vogt, he only mentioned it, 
briefly, in his textbook published in 1927 (Russell, Dugan & 
Stewart 1927) where full credit is given to Vogt. 
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scribes the physical properties of the tracks used in the 
analysis of the data, which are presented in Section §3. 
Section §4 deals with specific issues related to some inter¬ 
esting binary systems, and Sect §5 discusses the general 
results and provide new photometric [Fe/H] constraints. 
The conclusions are presented in Section §6. 

2. Evolutionary tracks 

Ideally tracks with a wide range of physical parameters 
(not only V and Z but also overshooting and convection 
parameters) should be used to properly test the different 
correlations between the parameters and the inferred ages 
and metallicities of the systems, so that proper statisti¬ 
cal uncertainties can be assessed. Previous analyses (e.g. 
P97, Allende Prieto & Lambert 1999, ROO, YOl) have re¬ 
lied on a given set, so cross-comparisons and systemat- 
ics are difficult to assess. It is unfortunate that there is 
presently no published set of tracks that fully fulfill this 
condition^ and we are limited to use different sets. In this 
study, we adopted three widely-used tracks computed by 

(1) the Geneva group (Schaller et al. 1992, Schaerer et 
al. 1993ab, Charbonnel et al. 1993, Mowlavi et al. 1998) 
which will be referred to collectively as the Geneva tracks, 

(2) the Padova group (Bressan et al. 1993, Fagotto et al. 
1994abc)^ and (3) the Granada group (Claret & Gimenez 
1992, hereafter CG92) evolutionary stellar models for a 
smaller range in metallicity than the Geneva and Padova 
ones. Table ^ summarizes the main physical ingredients 
that characterize these different tracks. We wish to em¬ 
phasize that our main concern here is to present and val¬ 
idate a general method whichever the models used, but 
not to test all the available (or the most recent) tracks: we 
only choose three different sets in order to understand and 
quantify possible systematic effects introduced by a given 
set. Kovaleva (2001) did a similar analysis using a subset 
of 43 systems and the Geneva and Padova tracks. We also 
describe in Table the older set from Hejlesen (1980ab, 
H80) because these tracks were often used previously for 
some of the systems studied here. They have a large mix¬ 
ing length theory (MLT) parameter (a= 2.0), do not in¬ 
clude overshooting or mass loss, and use the outdated Cox 
& Stewart (1970ab) opacities. For all these reasons we will 
only mention results derived from the H80 models as an 
historical guideline, because there is no reason to probe 
their validity (see comments on H80 related problems in 
e.g. Andersen et al. 1984). Some important points to note 
are the following : 

^ Grids of stellar evolutionary models for various values of 
initial helium abundances are available: Claret 1995, Claret & 
Gimenez 1995, Claret 1997 and Claret & Gimenez 1998 (re¬ 
ferred to as CG models) for metallicities Z=0.004, 0.01, 0.02 
and 0.03. However, overshooting and convection parameters 
are fixed in these models. 

^ We do not use the more recent Padova tracks (Girardi et 
al., 2000) because they do not include models for stars more 
massive than 7 M© and Z larger than 0.03, both relevant here. 


— Masses. For the 6 metallicities available in the Geneva 
tracks, the mass range is from 0.8 to 120 Mq (0.8 to 
60 Mq for Z = 0.1), while for the 6 metallicities in 
the Padova tracks the mass range is 0.6 to 120 Mq 
(restricted to the range 0.6 to 9 Mq for Z = 0.1). The 
CG92 grids cover masses from 1.0 to 40 Mq with 3 
different metallicities. The CG92 upper mass limit is 
large enough for the purpose of this work because the 
most massive component is a 32 Mq star (namely DH 
Cep). 

— Mass loss. The Geneva tracks use the mass loss rates 
given by De Jager et al. (1988) and Nieuwenhuijzen & 
De Jager (1990) for Pop. I stars. The Padova tracks 
set a mass loss rate only for stars more massive than 
12 Mq (from De Jager et al. 1988) with a metallicity 
dependence given by the prescription of Kudritzki et 
al. (1989). No mass loss by stellar winds is taken into 
account for the Z = 0.1 Padova tracks. Mass loss is also 
included in the CG92 tracks following Nieuwenhuijzen 
& De Jager (1990) and Reimers (1977) for cool giant 
stars. 

— Nuclear rates. The main difference between the Geneva 
and Padova tracks is the ^^(7(a,7)^®0 reaction which 
occurs during the helium burning phase, causing dif¬ 
ferences in the blue loop extension in the HR diagram. 
The Padova tracks include the nuclear reaction rates of 
Caughlan & Fowler (1988) which are somewhat smaller 
than the rates used by the Geneva group (Caughlan et 
al. 1985). 

— Overshooting. Although the nominal overshooting pa¬ 
rameter Ac of the Padova tracks is 0.5, this is just 
due to a different definition, and, as Fagotto et al. 
(1994a) explain, this value is equivalent to an over¬ 
shooting distance of Ac(Geneva) = 0.2 pressure scale 
heights. Claret & Gimenez (1992) have also adopted 
Ac = 0.2. To test the influence of overshooting on the 
results we also use Geneva tracks without overshoot¬ 
ing. However these models cover only the mass range 
[0.8-1.25] Mq and hence limit the application of the 
test to the few systems in this mass range. 

2.1. Ages 

One important definition which is not only semantic is the 
age of a star. Although some prescriptions use the start of 
the collapse of the parent cloud, we adopt here the con¬ 
vention of assigning a ‘zero age’ to the main sequence. 
The reason for this is twofold. First, the collapse phase 
is still poorly understood, and accretion processes before 
or during the Hayashi phase could well modify the time 
of arrival on the main sequence. Second, it seems simpler 
to follow the traditional definition of a ZAMS, a ’Zero- 
Age’ Main Sequence precisely when the radius of the star 
reaches a minimum value^. To avoid any possible confu- 

^ The usual prescription based on the H burning does not ap¬ 
ply for instance to massive stars, see e.g. Bernasconi & Maeder 
(1996). 
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Table 1. Basic properties of the Geneva, Padova, Granada and H80 tracks 




Geneva 



Padova 



Granada 


H80 “ 


Y 

Z 

Ref. 

Y 

Z 

Ref. 

Y 

Z 

Ref. 

Y 

Z 





0.230 

0.0004 

P94a 




0.2996 

0.0004 


0.243 

0.0010 

G92 










0.252 

0.0040 

G93b 

0.240 

0.0040 

P94b 




0.296 

0.0040 











0.196 

0.0040 


0.264 

0.0080 

G93a 

0.250 

0.0080 

P94b 













0.267 

0.0100 

CG92 

0.290 

0.0100 

Chemical 

0.300 

0.0200 

G92 

0.280 

0.0200 

P93 

0.280 

0.0200 

CG92 

0.380 

0.0200 

composition 










0.280 

0.0200 











0.180 

0.0200 








0.321 

0.0300 

CG92 

0.270 

0.0300 


0.340 

0.0400 

G93c 







0.360 

0.0400 











0.260 

0.0400 





0.352 

0.0500 

P94a 







0.480 

0.1000 

G98 

0.475 

0.1000 

P94c 







Opacities 

RI92, K91, IRW92, 
IR96, AF94 ^ 

IRW92, H77, CS70 " 

RI92, A92 (=AF94), 

H77, W90 

CS70 

MLT a 

1.6 

1.63 [1.5 for Z=0.1] 

1.5 

2.0 

A, 

0.2 

0.2 

0.2 

none 

Mass loss 

yes 

yes 

yes 

none 


“ Chemical composition {Y, Z) from Hejlesen (1980b). 

^ For Z = 0.001, 0.008 and 0.020, opacities are from Rogers & Iglesias 1992 (RI92). At low temperatures, from 6000 K to 2100 
K, they have completed the tables of RI92 with the atomic and molecular opacities of Kurucz 1991 (K91). For Z = 0.004 and 
0.040, the OPAL radiative opacities (Iglesias et al. 1992, IRW92) are used. These tables are also completed at low 
temperatures with K91. For Z = 0.1, opacities are from Iglesias & Rogers 1996 (IR96), along with the opacities from 
Alexander & Ferguson 1994 (AF94) at low temperatures. 

Because the OPAL tables of Iglesias et al. 1992 (IRW92) do not extend below 6000 K and above 10® K, the opacities at 
lower and higher temperature, respectively, are taken from the LAOL tabulations by Huebner et al. 1977 (H77) and Cox & 
Stewart 1970a,b (CS70). In general, the OPAL are significantly higher than the LAOL in two temperature ranges (a few 10® K 
and ~ 10® K). 

Opacities are taken from RI92, but for Teg < 6000 K or Tgg > 10® K they have used Alexander 1992 (A92, published later 
as Alexander & Ferguson 1994, AF94) opacities and Los Alamos Opacity Library - LAOL - (Huebner et al. 1977 [H77], Weiss 
et al. 1990 [W90]) respectively. 


sion, we hence define an ‘age’ r as the time elapsed after 
reaching the ZAMS, in logarithmic units : 

r = log(t - tzAMs) (1) 

where t is the actual (unknown) age and tzAMS the age 
when the star reaches the ZAMS (also unknown in abso¬ 
lute terms), both measured in years. We stress again that 
the actual age is meaningless unless a well defined ref¬ 
erence epoch is agreed on, and all the physical processes 
along the pre-MS evolution are fully understood and quan¬ 
tified. 

The Padova and CG92 tracks do follow the conven¬ 
tion of adopting tzAMS = 0) while the Geneva tracks do 
not. The dependence with mass and metallicity of the ra¬ 
tio tzAMS/txAMS (where TAMS is the Terminal Age Main 
Sequence i.e. central hydrogen exhaustion) for these tracks 
shows that for all practical purposes the non zero tzAMS 
constitutes a small correction of about 2% for evolved 


systems, but that is important for the youngest systems. 
All the ages derived in the present work are post-ZAMS, 
i.e. we explicitely substract tzAMS while dealing with the 
Geneva tracks. The H80 tracks seem to set tzAMS = 0 for 
all masses. 

2.2. Metallicities 

The set of tracks used in this paper do not have helium 
abundances Y independent of the metallicity Z (except 
those from H80), but rather follow the usual prescrip¬ 
tion of increasing Y for larger values of Z, roughly like 
Y = 0.23-1- (3.0± 1.0)Z (see e.g., Timmes et al. 1995). As 
shown in Table |^, the Geneva, Padova and CG92 tracks 
adopt different enrichment prescriptions. This means that 
the metallicity Z derived from a given set cannot strictly 
be compared with the Z derived from another set, since 
the different helium abundance will effectively mimic a 
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slightly different Z. This effect has to be kept in mind 
when comparing the derived metallicities from the hts. 

Usually, a simple relation is adopted between Z 
and [Fe/H]: [Fe/H]= \og{Z/X)—\og{ZQ/X q) (where we 
adopt solar abundance from Grevesse, 1997, priv. comm.: 
^0=0.017 and Xq=0.713). However for low metahicity 
stars, the abundance of heavy elements cannot be simply 
derived from spectroscopic measurements of [Fe/H] be¬ 
cause Oxygen (O) and a elements (Ca, Ne, Mg, S, Si) are 
overabundant in comparison to Iron (Fe). Therefore, be¬ 
cause we may have to consider low metallicities for some 
stars in our sample, throughout this paper we consider the 
correction proposed by Pagel (1996), 

^corrected = ^ X [0.638/a + 0.362] (2) 

where fa is a correction factor for O and the a elements: 

r 10°'4 for [Fe/H] < -1 

/„ = } io0.4x[Fe/H]= ^ [Fe/B] < -bO (3) 

[ 1 for [Fe/H] > +0 

3. Methodology 

3.1. The sample 

We have compiled a comprehensive list of 60 double- 
lined eclipsing binaries (EBs) which have the most ac¬ 
curate masses, radii and temperatures found in the liter¬ 
ature. The core of the sample is the large catalogue from 
Andersen (1991, hereafter A91), which is supplemented 
with 11 systems^. We also add three other systems with 
slightly lower accuracy criteria^ which can be used for a 
posteriori testing. A few systems have more recent and ac¬ 
curate measures than those listed by A91. We have used 
the updated parameters for 5 systems and for the EW Ori 
secondary component]^. It is worth noting that we kept all 
the systems from the Andersen list even EK Cep whose 
secondary component is in a pre-MS phase (which is not 
included in the stellar tracks that we used) because our 
method presented in j |3.3| can provide interesting results 
from the primary alone. 

One possible concern is that these systems may have 
undergone mass transfer episodes in the past. Although 
in some special cases strong constraints can be placed 
on the evolution of the system (e.g. Daems, Waelkens & 

® RT And (Popper 1994), AD Boo (Lacy 1997b), SW CMa 
(Lacy 1997c), AH Cep (Holmgren et al. 1990), CG Cyg (Popper 
1994), Y Cyg (Simon, Sturm & Fielder 1994), CM Dra (Viti 
et al. 1997), AG Per (Gimenez & Clansen 1994), V505 Per 
(Marschall et al. 1997), V3903 Sgr (Vaz et al. 1993, Vaz et al. 
1997), V526 Sgr and YY Sgr (Lacy 1997a). 

® AR Aur (Nordstrom & Johansen 1994a, mass accuracy 
about 4%), V477 Cyg (Gimenez & Quintana 1992, mass accu¬ 
racy: Myi~ 7% and M_b~ 4%), and DH Cep (Hilditch et al. 
1996, mass accuracy about 5%). 

^ V539 Ara (Clausen 1996), /3 Aur (Nordstrom & Johansen 
1994b), GG Lup (Andersen et al. 1993), EW Ori B (Popper 
1997), Al Phe (Milone et al. 1992), and DM Vir (Latham et 
al. 1996). 



Fig. 1. HR diagram for the sample of 58 binary systems 
with components more massive than 0.6 Mq. ZAMS se¬ 
quences computed by different authors for Z =0.02 are 
also indicated for comparison. 


Mayor 1997), in general this confuses the issue. However 
in our sample this is unlikely to be the case : (1) We 
have computed the Roche radii of the components and 
checked that they were much larger than the stellar radii 
(Lastennet 1998, from the formula of Eggleton 1983), and 
(2) most of the systems are still between the ZAMS and 
the TAMS, and their radii should not have evolved signif¬ 
icantly. The hypothesis of independent evolution without 
phases of mass transfer seems acceptable for every system 
in the sample, except perhaps for DH Cep. This system is 
an ellipsoidal variable, and the lack of eclipses makes any 
determination of the inclination very uncertain. Although 
Hilditch et al. (1996) and Penny et al. (1997) find that 
the system is detached, Burkholder et al. (1997) argue 
that the radii of the components are, according to their 
solution, very close to the Roche values. The uncertainties 
in the inclination of the orbit makes this system much less 
constraining. 

Figure shows the distribution of the sample on the 
HR diagram, where it is readily apparent that the systems 
are mostly sampling the evolutionary phases close to the 
main sequence. Also indicated are the ZAMS computed 
by different authors sX Z = 0.02, and the different pre¬ 
scriptions for overshooting and mass loss result in slightly 
different sequences for massive stars. Note that these se¬ 
quences are very different (by about 0.05 dex in Teff and 
0.5 dex in LjL q in the range of log Teff from 4.0 to about 
4.6) from the widely used sequence tabulated by Schmidt- 
Kaler (1982), which is no longer appropriate (Lastennet 
1998). 

Hereafter, each binary is associated with a number be¬ 
tween square brackets (e.g. /3 Aur [27]) corresponding to 
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its entry in all the tables and figures throughout this pa¬ 
per. 


3.2. Effective temperatures and luminosities 


The individual stars in the sample have typical relative 
errors smaller than 5% for the effective temperature^ and 
relative errors smaller than 10% for the luminosity. The 
interest of the luminosity is that it follows directly from 
the temperature and radius (Stefan’s law), and hence is 
independent of distance. Since analyses of these systems 
combine spectroscopy and photometry, the TeffS are likely 
to be more reliable than those in most of studies of single 
field stars. 

For these reasons, these 2 parameters were used di¬ 
rectly from the compilation of Andersen (1991) by dif¬ 
ferent authors (e.g. Lastennet, Lejeune & Valls-Gabaud 
1996, or Pols et al. 1997). However, because TeffS (and 
therefore luminosities) are far more indirect and inhomo¬ 
geneous quantities among the fundamental parameters, 
some efforts have been done to revise their determina¬ 
tions from photometric methods (e.g. Jordi et al. 1997, 
Lastennet et al. 1999a, Ribas et al. 2000, and Lastennet, 
Cuisinier & Lejeune 2002). 

We decided however to keep the original Tea values 


from the references given in §3.1 for simplicity because a 


comparison of the TeffS used by Ribas et al. (2000) and 
the ones used in this paper show a very close agreement 
for the 41 systems in common (~66% of our sample), the 
great majority (60 out of 82 stars) showing a difference 
below 2%. The Tea difference is less than 4% for all the 
components, except for 11 stars: both components of PV 
Gas [32], RS Cha [14], Al Hya [20], QX Gar [42] and IQ 
Per [33], and GZ GMa A [23]. Since our results are not 
reliable for the first two (bad fit for [32] and pre-MS stars 
for [14]) the Tea disagreement has no influence on our con¬ 
clusions, however the results derived for the latter systems 
([20], [23], [33] and [42]) have to be taken with caution be¬ 
cause the discrepancy in Tea would affect the luminosity 
as well, giving slightly different solutions. As suggested 
by Lastennet et al. (1999a, hereafter L99a), homogeneous 
determinations are highly preferred and to pick up some 
Teas in ROO and others in A91 and L99a would probably 
confuse even more the issue. Therefore for the purpose of 
this work and before a reliable homogeneous Tea determi¬ 
nation for all the systems of the sample, we decide to keep 


the original Teas (see §3.1). 


3.3. Fitting isochrones to binary systems 

There has been a variety of methods dealing with the 
problem of fitting isochrones to a given system. Besides 
the early attempts by Lastennet, Lejeune & Valls-Gabaud 
(1996), Pols et al. (1997) assumed that both components 

® Except in 6 less constraining systems where (j(Tea)/Teff is 
between 5 and 9.6%: V1031 Ori [28], V451 Oph [31], GG Lup 
B [36], EM Gar [45], DH Cep [53] and SW GMa [58]. 


are coeval and of the same metallicity and proceeded to 
minimize a functional which takes the observed masses, 
radii and temperatures and infers the best fit masses, 
common age and metallicity. Ribas et al. (2000) compute 
isochrones for each component, fixing each mass and grav¬ 
ity to the observed values, and proceed to compute a 
which minimizes the weighted difference between observed 
and predicted temperatures, as well as the difference be¬ 
tween the predicted ages of each component. Because the 
masses are fixed to the observed values, the x^ depends 
only on two remaining parameters, metallicity and helium 
content. In a recent paper. Young et al. (2001) use yet an¬ 
other formulation. First, they fix the metallicity to the 
solar value, and run specific evolutionary tracks of stars 
whose masses are given by the observations. Their x'^ func¬ 
tional minimizes the (sum of the squares of) residuals of 
luminosities and radii. Goevality is assumed, so that the 

depends only on age. 

A common property of these formulations is that it is 
not clear whether the observed quantities are used to their 
best discriminant power. Also the intrinsic correlations be¬ 
tween these quantities are not taken into account, and yet 
the very method of deriving these values, by relying on 
both photometric light curves and spectroscopic analyses, 
produces correlations among them. Another problem is 
that, for instance, masses are fixed (except in P97) to the 
mean observed values, so that even the observed disper¬ 
sions in mass are not taken into account. Obviously the 
dispersion in mass should have some effect on the final fits, 
but this is not considered (except in the P97 formulation). 
Although one could argue that using M, R and TeS (as 
P97 do) is the best way to proceed (no derived -and hence 
less precise- quantities, such as luminosity or gravity, are 
used), P97 must assume coevality and same chemical com¬ 
position. Stellar evolution theory is not therefore used to 
predict a priori values for each component. There cannot 
be possible discrepancies (which may well be the case, but 
it cannot be tested by this method). Likewise, ROO mini¬ 
mize the age difference between components, and derive a 
global metallicity and helium abundance for the systems. 

In a first formulation, we want then to test whether 
ages and metallicities can de derived for eaeh component 
independently. As a first step, we shall use here only the 
temperatures and luminosities of each of the components 
of the system. As it will be shown later on (see section 
§5.1) this does not make a difference in the final results. 
Note that because luminosity and temperature are not 
derived independently, the error bars in the HR diagram 
are highly correlated and strictly the error ellipse is always 
rotated (see e.g. Casey et al. 1998, for a brief discussion). 
To test the hypothesis that each component i {i = A, B) 
of a given system has the same age and metallicity, we 
derive them by minimizing a functional defined by 


X^{M„ti,Z,) = 


log L{track) — log Li 
cr(log Li) 

/ log Teff (track) - log Tern \' 
V Cr(l0gTeffQ ) 


( 4 ) 
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where the a are the 68% confidence level dispersions. 
The 3 sets of evolutionary tracks described in §2 are 
used to compute the isochrones by interpolating bilinearly 
at equal evolutionary phases in metallicity and age, and 
hence derive the values of L{track) and Teg(track). 

Because there are 3 parameters and only 2 data values, 
the problem is degenerate and one could find an infinity of 
solutions that fit the data. Among them we seek the min¬ 
imal values using the MINUIT package from the CERNLIB 
library (James 1994) and then explore the contours on 
a grid centred on the minimum minimorum. 

It turns out that, in practice, the mass is extremely 
well constrained and agrees with the measured value (see 
§ [I.4| ). This implies that a nearly unique solution (ti,Zi) 
ought to be found. In practice a whole set of solutions is 
found, which gives the empirical limitation to Vogt’s theo¬ 
rem. Let us illustrate the procedure taking as an example 
the Hya [34] system. Figure || shows the best fitting 
isochrones -using the Geneva tracks- for each component 
: they are not identical. The extent to which they are sta¬ 
tistically compatible depends on the confidence intervals 
for the metallicity and age of each of the components. The 
lower panels present the contours corresponding to the 
I, 2 and Sa bounds in this reduced distance, that is, the 
Xmm+ 1; ^ contours respectively. The degeneracy 

of the solutions is well illustrated here, as lower metallic¬ 
ity tracks fit a given component just as well at a slightly 
larger age, or a higher metallicity at a younger age. Note 
that at Ict, a factor of 2 in metallicity is allowed, whereas 
the uncertainty in the age is about 0.05 dex for compo¬ 
nent A and 0.3 for component B. The difference stems of 
course from the fact that A is more evolved, and hence 
better constrained than B. 

Note also that the 1 and 2cr contours on Fig. || are 
closed : within them the Vogt theorem is clearly violated, 
but since their areas are relatively small, one could argue 
that, in practice, the theorem is valid since a measure of 
the metallicity of one of the components would be enough 
to constrain very accurately its age. Nevertheless, how dif¬ 
ferent are isochrones distant la from the best solution ? 
Two examples are shown in Figure]^, where it is clear that 
even though the best fitting values for A and B are not 
the same (Fig. ^), at the Icr level they are statistically 
compatible. This is also apparent from the fact that the 
la contours of both components overlap : the ages and the 
metallicities of each component appear to be statistically 
compatible with each other. 

It remains to be seen whether the ages and metallici¬ 
ties of each component are the same as the ones inferred 
from the fitting of a single isochrone to the combined A-l-B 
system. Fitting the best isochrone to the combined system 
is done in a similar way, but in this case the x^ functional 
is obviously 


X (f sys : ^sys ) 


B 

E 

i=A 


log L{track) — log Li 
cr(log Li) 


logTesitrack) - logTegi V 
Cr (log Teffi) J 


(5) 


In this case the masses are no longer parameters (they are 
left free to take any value, just like the radii) and with 
4 data values and 2 free parameters, we expect to find 
a distribution with 2 degrees of freedom Using again 
MINUIT to find the true minimum value, we form the 
grid in the metallicity-age plane and compute the 1, 2 and 
3 ct bounds. Figure ^ shows a few isochrones within the 
la contour, and they are all compatible, at the Icr level, 
with the solutions found for each component (Fig. ^). Also 
illustrated are some isochrones chosen to show the effect 
of varying age, for a fixed Z, and changing metallicity, for 
a given age. In the case of x^ Hya, metallicity effects are 
clearly more important, and therefore a measure of their 
metallicities would constrain extremely well their common 
age within 0.05 dex (that is, about 11% at Icr). Note also 
that the solution found by Clausen and Nordstrom (1978) 
using the Hejlesen (1980) tracks is close to the 3cr contour 
(Fig. ^) and is therefore highly unlikely, as expected with 
models using the outdated Cox & Stewart (1970) opacities 

(cf. §i. 



Fig. 3. X^ Hya [34], individual components, Geneva 
tracks. For each component, the best fitting isochrones 
(from Fig. ^) along with the isochrones corresponding to 
the Icr level are indicated (dashed line : triangle; dot- 
dashed line : diamond, see Fig. ^). 

The irregularly-shaped contours are due to the differ¬ 
ent evolutionary speeds and short lived phases, and may 
complicate the interpretation in some cases. It is impor¬ 
tant to note, also, that since the contours are somewhat 
diagonal, the degeneracy between age and metallicity is 
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Fig. 2. Upper left panel: Zoom on the classical log Teff-log(L/Lo) HR diagram around the components of Hya [34] 
(Andersen 1991). In this example, an isochrone from the Geneva tracks is fitted to each component independently, 
without requiring equal ages or metallicities. The best fit isochrones are shown for each component (solid line : 
component A, dashed : component B). The isochrones derived from the HR diagram are also reported on the mass- 
radius diagram {upper right panel) and they show a very good fit to the very accurate masses and radii, even though 
these quantities were not used for the fit. Lower panels: Confidence regions in the (metallicity, age) plane at Icr (solid 
line), 2(7 (dashed line) and 3a (dot-dashed line) confidence levels for the primary (left panel) and the secondary (right 
panel) components. The triangles and diamonds indicate the {Z,logt) values used in Fig. 


lifted as soon as a measure of the metallicity is made. The 
contours for the system are much smaller than the con¬ 
tours for the individual components, as could be expected 
from the addition of information. 

The results of the fitting of, on one side, each compo¬ 
nent of the system, and, on the other hand, of the system 
on a single isochrone, are given in Table for all the sets 
of tracks which we have tested here. 


3.4. Masses and radii: the most stringent tests of the 
method 

To what extent these fitted isochrones are able to 
predict the observed values of the masses and radii? 
These are truly predictions from the theory, since the 
fit involved only luminosities and temperatures. Are the 


derived metallicities consistent with the observations? 
Metallicities will be discussed for the few cases where 
constraints exist in section §5.2, and here we analyse the 
predictions made by the fits for the entire sample to the 
very accurate masses and radii available (typical errors 
are smaller than 2% for mass and radius, except for 3 
systems where typical mass errors are between 4 and 7%). 


As clearly shown on Fig. and ^ the correct mass 
an radius are predicted for each individual star from 
isochrone fitting in the (log Tea, log L/Lq) HR dia¬ 
gram. The relative error for the derived masses of each 
component is given in Fig. for the Geneva tracks. 
The other tracks yield virtually the same results. Note 
again that the inferred mass comes from the fitting of an 
isochrone to the position of that component in the (log 
Teff, log L/Lq) diagram, without considering any other 
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Fig. 4. Hya system [34], Geneva tracks. Upper panel: 
Confidence regions in the (metallicity, age) plane at la 
(solid line), 2a (dashed line) and 3cr (dot-dashed line) con¬ 
fidence levels for both components. For the sake of com¬ 
pleteness, the best fit obtained by Clausen & Nordstrom 
(1978) using the H80 models is also shown (square with 
error bars). Lower panel: Corresponding isochrones on 
the HR diagram (error bars from Andersen 1991) as in 
Fig.|^. The best fit isochrone obtained with the Geneva 
tracks and defined by {\ogt,Z) = (8.275, 0.014) is indi¬ 
cated by the solid line. The 4 points on the la contour 
above correspond to the 4 isochrones showed in the lower 
panel. 

constraint. A linear regression yields AM/Mmeasured = 
(Afderived Af^^^easured)/Afmeasured — ( 4.2 zb 5.7) 10 
(6.35±1478.3) 10“® Mmeasured, that is, both slope and sys¬ 
tematic trend consistent with zero. The average relative 
error in mass is —0.54%, with a dispersion of 1.7%. This is 
a remarkable achievement for the theoretical tracks, since 
over this mass range we are probing stars which are almost 
fully convective to fully radiative stars. 

Only a few components (4 out of 116) present a signif¬ 
icant discrepancy between the predicted mass and the ob¬ 
served one, given the error in the observed value, in Fig. 
These are component B of AD Boo [57], both components 
of PV Cas [32] and IQ Per [33]. AD Boo B yields in fact 
a bad fit, and one cannot expect to infer a correct mass 
or radius from a bad fit. However, a more recent determi¬ 
nation of the secondary component mass (Mb = 1.20±0.03 
Mq, Popper 1998) than the one we used (1.237±0.013 Mq, 
Lacy 1997b) gives a better agreement with our theoretical 
predictions from the Padova (1.188±0.024 Mq) or Geneva 
models (1.155^0 022 excellent agreement 

with the Granada tracks: 1.205±0.025 Mq. 


Fig. 5. Relative error in mass, AM/Mmeasured = 
(Afderived Afiaeasured)/Afmeasured 5 inferred from the fitting 
of Geneva isochrones to the position of each component 
on the HR diagram. The dotted line is AM = 0 and the 
dot-dashed lines are linear regressions to the data points. 





R(measured) [Rg] 

Fig. 6. Relative error in radius, Ai?/i?nieasured = 
(Rderived 77measured)/77measured: inferred from the fitting 

of Padova isochrones to the position of each component 
on the HR diagram. The dotted line is AT? = 0 and the 
dot-dashed lines are linear regressions to the data points. 

The problem is the same for PV Cas [32], which is 
badly fitted by the tracks we have used. Pols et al. (1997) 
also obtained bad fits, even at the upper limit in Z they 
explored [Z =0.03), and despite the fact that we explore a 
larger range (up to Z>0.04) we still do not manage to get 
a correct fit. Popper (1987) failed also to fit it with 3 other 
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sets of models: the predicted Tes is too hot by ~1700 K. 
The revised ROO Tegs, lower by about 700 K, would give 
a better agreement with solar metallicity models, but as 
suggested by Young et al. (2001) this system may be still 
in a pre-MS phase, which is not taken into account in the 
tracks studied here. 

The third, IQ Per [33] doesn’t produce a too bad fit, 
and, as we discussed in [ j3.2| , adopting the ROO Tegs would 
give a better agreement. 

As far as the radii are concerned. Fig. || shows 
the same trend as for the masses. A linear regres¬ 
sion yields AR/Rmeasured = ( — 2.2 ± 1.2) 10“^ -I- (0.8 ± 
0.4) 10“^i?nieasured- The slightly positive trend is produced 
by a few outliers which however are not separated by the 
measured values by more than la. Overall the mean rel¬ 
ative error is 0.03% with a dispersion of 3.8%, again a 
remarkable achievement for the theoretical tracks, since 
we are probing stars which are almost fully convective to 
fully radiative stars. 

Hence, the positions of these components in the HR 
diagram yield extremely precise and accurate predictions 
on their radii and masses, with dispersions of below 4% 
around a mean error below 0.5% both in radius and mass. 


3.5. Tests and limitations 

We should however be cautious before taking these re¬ 
sults at face, because one has to realize that these errors 
are only a combination of internal errors plus a fraction 
of external ones. Not all the systematic errors have been 
taken into account, nor the covariance terms in the inter¬ 
nal, measurement errors. We will briefly mention here four 
important effects that one has to keep in mind when us¬ 
ing these results: the Helium abundance, the availability 
of tracks with a wider range of Z, overshooting, and the 
improved accuracy of the measures. 

(1) One obvious test case of a single component is 
the nearest star, the Sun, for which the effective temper¬ 
ature is known to 0.04% and the luminosity to 0.26%. 
Figure ^ shows the results of fitting the Sun as a (sin¬ 
gle) component using the Padova and the Geneva tracks. 
The very signihcant difference, at first sight very surpris¬ 
ing, can mainly be ascribed to the larger He abundance 
used by the Geneva group in comparison to the Padova 
tracks as shown in Tab. (Charbonnel & Lebreton 1993). 
The effects of the equation of state and opacities (Yildiz 
& Kiziloglu 1997) or diffusion (Morel et al. 1997, Weiss & 
Schlattl 1998) are comparatively much smaller. In the con¬ 
text of our work, this difference is a simple caveat not to 
use the published tracks blindly, and also calls for the cal¬ 
culation of a large grid of models with uncorrelated vari¬ 
ations in the main physical parameters. The derived ages 
and metallicities are therefore track-dependent, and the 
values for (T(Z),cr(T) that we quote are necessarily lower 
limits because the systematic effects introduced by the 
tracks cannot be taken into account. The best we can do 
is to compare the results obtained with different tracks to 


assess the statistical robustness of the values derived with 
this procedure. 

(2) For systems where the components are quite close 
to each other in the HR diagram, the constraints coming 
from the fitting of the system are weaker since there is 
effectively only one star as a constraint, hence increasing 
the degeneracy. Note that even if the contours produced 
by each component are similar, the best htting values can 
be very different. Sometimes this difference arises from the 
fact that the best value is close to the boundary in metal¬ 
licity allowed by the tracks. In several cases (see Table 
the htting indicates large metallicities that can only be 
reached with the Padova and Geneva tracks. The best hts 
obtained are not necessarily good hts. 

(3) The 3 sets of tracks used here have the same pre¬ 
scription for core overshooting (see Tab. |^). Would tracks 
without overshooting produce worse or better hts? For in¬ 
stance, P97 argue that a few systems can only be htted 
with models which include overshooting. We have used 
Geneva tracks without overshooting for this test, and as 
shown in Tab. we hnd no signihcant differences be¬ 
tween the Geneva models with and without overshooting. 
However, the power of the test is limited by the small mass 
range (and hence the number of systems) in common for 
these two sets of models, and by the fact that most of these 
systems yield bad hts (see §0). We cannot conclude, on 
the basis of these tracks and this small sample, whether 
overshooting is preferred or not. This is yet another source 
of systematic uncertainty that is difficult to quantify. 

(4) To what extent the age-metallicity degeneracy can 
be partially lifted by increasing the accuracy of the mea¬ 
sure^ A good example is provided by V539 Ara [41] 
(§ 4.5.2 and Fig. where the measures from Clausen 
(1996) are a factor of about two better than the previ¬ 
ous ones from Andersen (1983). The degenerate areas in 
the {t, Z) plane are substantially reduced, although large 
metallicit y soluti ons are allowed. See also the case of /3 
Aur [27] ([ j4.4.10| and Fig. ^) for a similar example. Once 
again, the most contraining systems are those with evolved 
components, even though the precision in the parameters 
could be worse. 


4. Analysis of individual systems 

We have applied the techniques developed in §3 to the 
entire sample of 60 systems. The results are summarized 
in Table ^ for 58 systems (from [2] to [59]) using the 3 
different sets of tracks described in §2. The two remain¬ 
ing systems ([la] and [lb]) are discussed in §^. The col¬ 
umn marked ’T’ indicates the track used for the ht. Note 
also that all ages were corrected for the finite age of the 
ZAMS in the case of the Geneva tracks, following Eq. (1). 
Although the best ht values may appear discrepant, it is 
essential to use the minimum and maximum values tab¬ 
ulated (at the n-th conhdence level, as indicated) to test 
for the statistical signihcance of the differences. 

It is beyond the scope of this paper to analyse each 
binary system in detail, hence in the following we only 
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Fig. 7. Using the Sun as a test case for the Padova (diamond and open star) and Geneva (circle and triangle) tracks. 
The adopted temperature and luminosity are 5777 ± 2.5 K and (3.844 ± 0.01) x 10^^ erg s“^ (Weiss & Schlattl 1998). 
A bad fit (open cross) obtained with the Geneva models but with parameters close to the Z-\ogt Padova solutions is 
also shown for comparison. 


discuss a few systems which either can be compared with 
previous studies or present problems. Except pre-MS stars 
(§4.1) were good fits cannot be expected with post-MS 
models, bad fits were obtained in the HR diagram for stars 
less massive than 1.2 Mq, EW Ori [3] being the most 
massive one. They are mainly discussed in §4.3. 


4.1. Binaries with (at least) a pre-MS component: EK 
Cep [17] and RS Cha [14] 

As expected, none of the three sets of tracks are able to 
fit the system EK Gep [17], because the secondary is a 
pre-main-sequence stain and the tracks used do not in¬ 
clude these phases. Nevertheless, since our method can use 
information for individual stars, we can predict from the 
primary alone a metallicity, which turns out to be around 
solar. 

RS Cha [14] is also identified as a pre-MS (Mamajek et al. 

2000) : pre-MS tracks achieve a good fit (Palla & Stabler 

2001) . RS Cha will not be discussed any further due to 


® As suggested by Popper (1987), and confirmed by Martin 
& Rebolo (1993), and Claret, Gimenez & Martin (1995). 


the uncertainty on its Tea (see §3.2). 

Finally, as suggested by Nordstrom & Johansen (1994a, 
hereafter NJ94a), AR Aur B [48] may also be a pre- 
MS as well. Tab ^ predicts a metallicity slightly sub¬ 
solar (Z~0.014±0.003) for the system, in agreement with 
the NJ94a study. We derived from synthetic photometry 
(L99a) a lower limit on the metallicity of the primary 
(2’>0.017), only marginally consistent with the theoreti¬ 
cal results. This discrepancy may come from the chemical 
peculiarities of both components (Zverko et al. 1997, and 
references therein). 


4.2. Binaries with VLM star components (M<0.7 
Mq): YY Gem [la] and CM Dra [lb] 


Neither YY Gem [la] nor CM Dra [lb] provide a test for 
Geneva, Padova and Granada tracks, since their compo¬ 
nents are very low mass (VLM) stars, with masses below 
0.6 Mq. For these two systems only, we have used instead 
the evolutionary tracks from Baraffe et al. (1995) which 
are computed at logt = 10 yrs. Assuming the position of 
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Fig. 8. YY Gem [la] and CM Dra [lb]. Data points for 
CM Dra come from Viti et al. (1997). The circle gives the 
position of YY Gem from Andersen (1991, both compo¬ 
nents at the same position), to be compared with Leung 
& Schneider (1978, error bars for each component) and 
Torres & Ribas (2002, bold error bars, both components at 
the same position). The isochrones of Baraffe et al. (1995) 
and Brocato et al. (1998) are for a solar metallicity and an 
age of 10 Gyr. Two analytical ZAMS (2’=0.02 and 0.03) 
from Tout et al. (1996) are also shown for comparison. 

YY Gem from A91^ (see Fig. these tracks yield a best 
fit for YY Gem at ^ = 0.019. This is consistent with the 
solution found by Chabrier & Baraffe (1995), indicating 
also that YY Gem [la] could be anywhere between the 
upper limit of 10 Gyr and the end of the pre-MS contrac¬ 
tion phase. Since this upper limit may be too old (Torres 
& Ribas, 2002 estimated a younger age of 370 Myr), 
we refer the interested reader to the detailed analysis of 
YY Gem by Torres & Ribas (2002). We also note that us¬ 
ing the older Tefj-L parameters from Leung and Schneider 
(1978) (also shown on Fig. leads to virtually no con¬ 
straint due to the large error bar^. 

Torres & Ribas (2002) adopted Teff=3820±100 K (see 
Fig. ^ from several empirical calibrations. However, 
adopting the colour index values listed in their Tab. 6 
for V—Ic, V—K and (R—I)c, assuming solar abundances 

The comprehensive study of YY Gem by Torres & Ribas 
(2002) provides a new position in the HR diagram which how¬ 
ever leads to no significant differences. 

As shown by Torres & Ribas (2002), the apparently good 
fit in the HRD for YY Gem using the Baraffe models is de¬ 
ceiving because they fail to predict the correct radius of low 
mass stars for the measured mass, by 10% or more. At least 
in this particular case, the use of L and Tea gives the wrong 
impression. This illustrates the importance of the prescription 
proposed in this paper: 1) to fit the classical HRD and then 2) 
to check the predicted masses and radii. 


and the very accurate surface gravity of YY Gem from 
Torres & Ribas (2002), and neglecting the reddeninj^ 
the BaSeL models (Lejeune et al. 1997, 1998) indicate 
cooler TeffS (between 3605 and 3642 K to be compared 
with the coolest temperature that Torres & Ribas (2002) 
obtain: 3719 K from V—Ic). Since the BaSeL models pro¬ 
vide very good agreement with empirical calibrations in 
these colours (cf. Fig. 2 of Lejeune et al., 1998), we sus¬ 
pect that the TeffS of YY Gem still need further revisions, 
but the temperature issue for cool stars is still controver¬ 
sial to give a definitive Teff . 

In the case of CM Dra [lb], whose metallicity is compat¬ 
ible with the solar value (Gizis 1997) or lower (Viti et al. 
1997, see also Metcalfe et al. 1996), its position in the HR 
diagram is very well constrained (Viti et al. 1997) for such 
a VLM binary, and indicates a metallicity close to solar 
with the models of Baraffe et al. 1995 (see Fig. ||). When 
using other sets of tracks (from Brocato et al. 1998, or the 
analytical ZAMS of Tout et al. 1996), the results are sub¬ 
stantially different and a. Z = 0.02 composition is clearly 
inconsistent with the data (see Fig. ^). A larger sample is 
needed in these range of masses to confirm these tests, and 
the discovery of an M-dwarf eclipsing binary, GJ 2069A 
(Delfosse et al. 1999), with components of about 0.4 Mq 
( which places GJ 2069A between CM Dra and YY Gem), 
is a first step towards better observational constraints in 
the lower mass end of the main sequence (see Segransan et 
al. 2000 for a recent list of accurate masses of VLM stars). 

4.3. Systems with at least one component in the 

0.7-1.1 Mq mass range: RT And [47], HS Aur [2], 
CG Cyg [46], and FL Lyr [4] 

We obtain systematically worse fits for all these systemsf^ 
The bad fits however present two different aspects. While 
both components of HS Aur [2] can not be fitted (or at 
least only by very enriched models), the systems RT And 
[47], CG Cyg [46] and FL Lyr [4] have a secondary com¬ 
ponent far too cool to be matched by the same isochrone 
as the primary]^. This may reveal a different behaviour. 
Under the optimistic hypothesis that the location of both 
components of HS Aur in the HR diagram are perfectly 
reliable, Fig. ^suggests that the models are unable to cor¬ 
rectly predict its properties: either we have strong difficul¬ 
ties to fit both components of HS Aur (Geneva models) or 
we obtain very questionable metal-rich solutions (Padova 
models). For illustrative purposes, Fig. ^ shows that the 
Geneva tracks do not produce confidence areas at the Icr 
level, and most solutions are extremely far away (3 ct) from 


YY Gem is a nearby system (~16 pc) according to its 
Hipparcos parallaxe, so we assume a colour excess E(B—V)=0. 

Some of their components are not massive enough to be 
fitted with the Granada tracks, so only the Geneva and Padova 
tracks are used here. 

The same difficulty appears with the Gambridge group 
models by fitting simultaneously their effective temperatures, 
masses and radii (see Pols et al., 1997). 
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Fig. 9. HS Aur system [2]. Left panels : Confidence regions in the (metallicity, age) plane. Right panels : HR diagrams. 
The top panels use the Geneva tracks, while the lower ones use the Padova ones. The best solutions obtained are 
indicated by circles. Keeping the same metallicity as the best solution, we also show two other points at 3 (t from the 
best solution, with respectively a smaller age (diamond) and a larger age (open star). Another point (triangle) is the 
solution with the lowest metallicity allowed at a 3a confidence level. The Geneva isochrones clearly do not fit the HS 
Aur system, no matter what age or metallicity are used. In contrast, the Padova tracks allow a successful fit at about 
2.5 times the solar metallicity. 


the position of the system. In contrast, the Padova tracks 
succeed in fitting the system, but at a metallicity which is 
more than 2.5 times solar (Z~ 0.04-0.05). This anomaly 
was also found by Pols et al. (1997): they obtained high Z 
values for all systems with stars below 1 Mq . The Padova 
tracks indicate that HS Aur is at the same time quite old 
and metal rich (Fig. ||, lower panels): ~[Fe/H]=0.4-0.5. It 
remains to be seen whether a spectroscopic analysis con¬ 
firms this predicted large metallicity. 

Since the results of Pols et al. (1997) and the present study 
rely on the same Teff determination, a revision may well 
be needed (in particular the Tefjs of the secondary compo¬ 
nents) before definitive conclusions. Unfortunately, while 
L99a, ROO, and Lastennet, Guisinier & Lejeune (2002) 
carefully re-derive Tefjs for many EBs from photometric 
calibrations, none of these works study the stars in ques¬ 
tion because individual uvby photometry would be neces¬ 
sary. 

We now have to examine FL Lyr [4], GG Gyg [46] and 
RT And [47]. Besides a Tea revision, possible explana¬ 


tions of the disagreement may come from mass transfer 
and starspot activity because these 3 systems are known 
to be active stars. As already discussed in §3.1 none of 


these stars overflows its Roche lobe (see also Lastennet et 
al., 2002 for more details). Yet in RT And, the face-to- 
face position of the spots on the surface of both compo¬ 
nents may indicate the possibility of a mass transfer from 
the primary to the secondary component through a mag¬ 
netic bridge connecting both active regions (Pribulla et al. 
2000). If this is confirmed, the RT And components could 
not be assumed anymore to evolve like single stars, and 
we could not expect to reproduce their evolution with the 
models considered in this paper. Enhanced surface activ¬ 
ity is also present in FL Lyr and CG Cyg. These spots lead 
to quite distorted light curves and consequently less well- 
determined photometric elements, and hence may provide 
an explanation of the disagreement we have found. 

If the disagreements appear to remain even after fur¬ 
ther analyses, they will give additional examples to the 
dilemma pointed out by Popper (1997) who remarked 
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that for systems with masses between 0.7 and 1.1 Mq 
it was difficult to place both two components on the same 
isochrone. The same happens in a Mass-Radius diagram 
as shown by Clausen et al. (1999). 

Other binaries show similar problems in this mass rangePI, 
as briefly reviewed by Lastennet et al. (2002). However, 
the systematic differences between models and binaries of 
this mass range are not observed in CD Tau C, a solar 
mass companion of the triple system CD Tau. Ribas et al. 
(1999) obtained a perfect fit of the three components with 
a single isochrone, and the issue is still open to debate. 
New accurate data for low mass eclipsing binaries should 
help to work out the discrepancies (Clausen et al., 2001; 
see also Lastennet et ah, 2002 and references therein). 

4.4. Systems with both components in the 1.1-3. Mq 
mass range 

4.4.1. Al Phe [5] 

Al Phe is one of the most interesting systems in the sample 
since the coolest component is also away from the main- 
sequence, and observational constraints on Z are avail¬ 
able. A spectroscopic determination of metallicity from 
high resolution CCD spectra (Andersen et al. 1988) gives 
[Fe/H] = -0.14 ± 0.1, i.e. Z = 0.012 ± 0.003. Using 
tracks from VandenBerg (1983) they derived an age of 
4.1 ± 0.5 X 10® yr from isochrones with an helium abun¬ 
dance Y = 0.27 ± 0.02p|. Our 3 sets of tracks can only 
fit the system at slightly lower metallicities, as shown in 
Fig. I^for the case of the Geneva tracks. These results are 
consistent with P97 and ROO. The combination of the up¬ 
dated parameters for the system from Milone et al. (1992) 
with the heavy element abundance from Andersen et al. 
(1988) constrain extremely well the system, with uncer¬ 
tainties in age smaller that 0.010 dex at Itr and 0.014 
at 3tT. The best fit solutions are entirely compatible with 
both data sets. The results obtained with the Granada 
tracks {Z = O.OlliO.OOl, 9.64 < logt < 9.67 yr, all 
Icr) and with the Padova ones {Z = 0.008to'ooo 9.59 

< logt < 9.63, Icr) are both in good agreement with the 
ones deduced from Geneva and the metallicity range from 
Andersen et al. (1988), although the range allowed in age 
is a bit larger. 


The Padova isochrones do not match the Hyades eclipsing 
binary V818 Tauri in the mass-radius diagram (thus, with¬ 
out using any information on Tefi), and it appears difficult to 
match the mass of the secondary of the nearby visual binary 
85 Peg. 

This abundance is much smaller than the one found by 
VandenBerg & Hrivnak (1985), Y = 0.33-0.44, although for 
a similar age (t — 3-4 x 10® yr). VandenBerg & Hrivnak 
(1985) adopted [Fe/H] = 0.17 ± 0.20 (from the observed mi 
index, the Crawford (1975) calibration for F—type stars and 
the (5, [Fe/H]) relationship of Crawford & Perry (1976)), which 
is clearly far too large in comparison to the determination of 
Andersen et al. (1988). 



Fig. 10. Al Phe [5], Geneva tracks. The upper panels 
use the parameters from Andersen et al. (1988), while 
the lower ones take the more recent data by Milone et 
al. (1992). Circles (solid line isochrones) are the best 
solutions. In the upper panel the triangle (dashed line 
isochrone) is at Icr from the best solution. The spectro¬ 
scopic metallicity {Z =0.012±0.003) from Andersen et al. 
(1988) is indicated with vertical lines. 


4.4.2. UX Men [6] 

This system was studied by Andersen et al. (1989) who 
derived a best fit isochrone of age (2.7±0.3) x 10® yr and 
Y = 0.27±0.01 using the VandenBerg (1985) tracks fixed 
at Z = 0.019, the observed value. There are several mea¬ 
sures of the metallicity]^ Since they are consistent with 
each other (except that the Clausen & Grpnbech value 
was not corrected for reddening), we adopt in Fig. the 
more precise spectroscopic determination (Andersen et al. 
1989). 

The Geneva tracks (with overshooting) give two coeval 
components with the same metallicity^. Even though the 
ages are similar to the one derived by Andersen et al. 
(1989), the metallicity derived from the best fit is 50% 

(1) Clausen & Crpnbech (1976) find a small value, [Fe/H] 
= —0.15 from uvby/3 photometry; (2) Kobi & North (1990) ob¬ 
tain [Fe/H] = —0.04; (3) From their own Stromgren uvby/3 ob¬ 
servations and correcting for reddening, Andersen et al. (1989) 
deduce from Nissen’s (1981) calibrations a metallicity of [Fe/H] 
= —0.05±0.15; (4) Andersen et al. (1989) also obtained from 
high resolution CCD spectra [Fe/H] = 0.04±0.10 (i.e. Z — 
0.019±0.004 for an assumed Zq = 0.0169, VandenBerg 1985). 

The primary has {logt,Z)A = (9.42, 0.033) whereas the 
secondary has (logt, Z)s = (9.44, 0.030) (see Table ^ for the 
uncertainties). 
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Fig. 12. DM Vir system [8]. Confidence regions using the Latham et al. (1996) data and the Geneva {upper panels), 
Padova {middle panels) and Granada {lower panels) tracks. Best fit values (circles) are indicated in each panel, along 
with the corresponding isochrone (solid lines) in the HRD. Other isochrones, associated with the positions marked 
by triangles or diamonds on the left panels, are also shown in the HRDs. The fit obtained by Andersen, Clausen & 
Nordstrom (1984b) with the H80 tracks (open star) and by Latham et al. (1996) with the Granada tracks (square) 
are also indicated. 


larger than the spectroscopic/photometric valuep|. This 
is also supported by the fitting of the combined system 
(Fig. M). The Padova tracks give a very similar result, 
{\ogt,Z) = (9.257, 0.034), but the Granada models and 
P97 yield a slightly smaller metallicity ((logt, Z) = (9.405, 
0.027) and (9.365, 0.024) respectively). However, at the 
Icr confidence level a value as small as Z =0.022 (0.025) 
is allowed by the Geneva (Padova) models. We suspect 
that a slight increase in the Helium abundance may rec¬ 
oncile the range inferred from its position on the HR di¬ 
agram with the spectroscopically determined metallicity. 
Alternatively, a revision of the Teg may be another expla¬ 


nation. This will be discussed further in §5.2 


An even more metallic abundance is suggested by the ROO 
extrapolated solution: Z—0.039. 


4.4.3. DM Vir [8] 

Accurate data from Latham et al. (1996) have superseded 
the values listed by A91 for this system. Figure ^ shows 
the best fits obtained with the Geneva (upper panels), 
Padova (middle panels) and Granada tracks (lower pan¬ 
els), and illustrates the danger of using a single set of 
tracks to get all the possible solutions. The limited Z range 
in the Granada tracks lead Latham et al. (1996) to a best 
fit at the Z limit of the Granada models, but which is by 
no means unique, as shown in Fig. |^. A similar good fit 
was also found by Pols et al. (1997) at Z~0.03, which is 
the upper limit of the models they used. This system was 
not selected by ROO. 

Although the age is very well constrained by all sets of 
tracks (about O.ldex at Icr, see Fig. |l^, left panels), solu¬ 
tions with metallicities ranging from 0.025 to 0.06 are not 
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Fig. 11. UX Men system [6], Geneva tracks. Best fit (open 
circle), Icr points (triangles and diamonds), and a 2a 
point (open star) situated within the observed metallic- 
ity range. The vertical lines mark the observed value of Z 
at 0.019±0.004 (Andersen et al. 1989), while the best age 
(open square) by Andersen et al. (1989) was derived using 
the VandenBerg (1985) models computed for Z = 0.019. 
Results from Pols et al. (1997) (error bars at Z =0.024) 
and the extrapolated solution of Ribas et al. (2000) (ver¬ 
tical dotted line at Z =0.039) are also shown for compar¬ 
ison. The last one is unreliable because well beyond the 
upper Z-limit covered by their models (Z=0.03). 


statistically rejected at the Icr level^. Since toi(DM Vir) 
— TOi(Hyades) = 0.001 (Andersen, Clausen & Nordstrom 
1984b), this suggests that DM Vir has a metal content 
similar to the one in the Hyades. Latham et al. (1996) 
also take [Fe/H] = 0.12±0.12, although a detailed spectro¬ 
scopic analysis is needed. This would imply (with Yq = 
0.27) Z = 0.023±0.006, yet the best fits obtained here 
with the 3 sets of tracks favour much larger metallicities, 
especially the Padova tracks. Hence, if this metallicity is 
confirmed, the contours of Fig. show that the Padova 
tracks would fail to fit the system in the HR diagram, 
while the Geneva and Granada ones would be consistent, 
provided the He abundance is as high as the one used by 
these models (cf. Tab. p. 


This was also found by Andersen, Clausen & Nordstrom 
(1984b) using the H80 models : solutions as separated as 
(logt = 9.60, Z = 0.02, Y = 0.18) and (logt = 9.60, Z = 
0.04, Y = 0.26) being good fits, the solar metallicity solution 
implying however a rather unlikely low helium abundance. 


4.4.4. RZ Cha [10] 

The components of RZ Cha are known to be both evolved 
and older than 2 Gyr (Jprgensen & Gyldenkerne 1975). 
Andersen et al. (1975) found logt = 9.301 with a prelim¬ 
inary version of H80 tracks, while we obtain (logt, Z) = 
(9.32to;o3, O.OISIq;™^) with the Geneva tracks. As illus¬ 
trated in Fig. 1^, the contours resulting from the fit with 
the Padova and Granada models are also very similar^ 
Yet these solutions for Z are systematically smaller than 
the value of [Fe/H] = —0.02 ± 0. 150 which is almost solar 
metallicity, Z = 0.020/jg QQg. In fact, the Hyades metallic¬ 
ity adopted by Jprgensen & Gyldenkerne was an over¬ 
estimation. More recent determinations (e.g. Cayrel de 
Strobel et al. 1997, Perryman et al. 1998) cluster around 
[Fe/Hjnyades = 0.14 ± 0.05, so that we derive for RZ Cha : 
[Fe/H] = -0.08 ± 0.05 or Z = 0.014 ± 0.002. This metal¬ 
licity is in very good agreement with the best fits obtained 
by the 3 sets of tracks, the Geneva tracks being slightly 
better (see Fig. |^). The extremely good constraint on 
the age (less than 0.04 dex at la) arises from its evolved 
position on the HR diagram. Despite the fact that both 
components of RZ Cha are on the same point in the HR 
diagram , hence a less than optimal case (Fig. |^), it is 
precisely because its position lies by the TAMS, that one 
can derive good constraints on both metallicity and age. 

4.4.5. PV Pup [12] 

The proximity of the two non-evolved components of this 
system in the HR diagram does not allow to define a very 
accurate solution: the closer the two points on the HR di¬ 
agram, the larger the contours. The Geneva tracks give a 
best fit at (logt, Z) = (7.53, 0.040), but metallicities as 
low as Z = 0.015 and ages as large as logt = 9.1 are not 
excluded at the 3a confidence level (Fig. j^. As indicated 
by Vaz & Andersen (1984), the H80 models allow good fits 
at (logt, Z, Y) = (9., 0.02, 0.18) or at (8.60, 0.04, 0.26), 
yet another example of the age—metallicity degeneracy. 
Nissen’s (1981) calibration of the Stromgren j3 and <5mi 
indexes gives Z = 0.017 (this value is a vertical line in the 
iso-y^ diagrams in Figs. [^ and For clarity. Fig. 14 
shows the results with the Geneva models and Fig. 15 
with Padova models. They both indicate that better fits 
are obtained with larger metallicities (in agreement with 
Pols et al. 1997: Z >0.028), yet isochrones at Z = 0.017 
are only 2a away from the best fits. The revised ROO Tefjs 
are the same ones that we use, so our high-metallicity re¬ 
sults are unchanged. Alternatively, the simultaneous (Teff, 
[Fe/H]) solutions derived from Stromgren photometry for 

(logt,Z) = (g.SOtg g^, 0.01812 ggl) and (logt,Z) = 
(9.3412^071 0-0f’7lto.oo4) Padova and Granada models re¬ 
spectively. These results are in a very good agreement with 
P97 (their overshooting models). For the sake of completeness, 
we note that this system was not studied by ROO. 

22 Derived from Ami = mi(RZ Gha) — mi (Hyades) = 
0.017 and choosing [Fe/H]Hyades = 0.20 ± 0.15, Jprgensen & 
Gyldenkerne (1975). 
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Fig. 13. RZ Cha system [10]. The upper panel superposes 
the la contours from the 3 sets of tracks, along with the 
best fits from Geneva (open circle), Padova (square) and 
Granada (triangle) tracks. The metallicity range found by 
Jprgensen & Gyldenkerne (1975), corrected for the up¬ 
dated Hyades metallicity, is indicated by the vertical lines. 
The best fit isochrones from the three sets of tracks are 
shown in the lower panel. 

both components of PV Pup (Lastennet et al. 1999a) give 
weak constraints on the metallicity (but an upper limit 
of [Fe/H]<0.24, i.e. Z <0.029, at 2-tT) and a rather well- 
defined range of Tea (log Tea between 3.86 and 3.88). This 
would shift the position of PV Pup in the HR diagram to¬ 
wards larger values of Tea (and hence L) and consequently 
would exclude all the high-Z solutions. 

4.4.6. TZ For [18] 

TZ For is a rare occurrence of a system with a sub—gi¬ 
ant and a giant, and hence potentially one of the most 
constraining systems. Moreover, there are spectroscopic 
data for this binary: Andersen et al. (1991) give [Fe/H] 
= -|-0.1±0.15^ The primary of TZ For is too evolved 
(core helium burning phase, as suggested by Claret & 
Gimenez, 1995a and Pols et al., 1997) to be matched by 
Granada models, hence only results derived for the B com¬ 
ponent are given in Tab. |^. The agreement around solar 
metallicity is quite satisfactory for the Geneva models, 
but seems difficult to match with the Padova tracks: the 
lower limit of the Padova models at the Icr level is only 
marginally consistent with a solar metallicity. In spite of 
the better agreement obtained with the Geneva models, 

i.e. Z = 0.021^0 QQg (with Zq = 0.0169, VandenBerg 1985) 
or Z = 0.0231 q;qJ° (for Zq = 0.0189, Maeder & Meynet 1988). 


Fig. 14. PV Pup [12], Geneva tracks. Best ht (open circle 
at (logt,Z) = (7.53, 0.04)), Icr points (triangle at Z = 
0.04, and diamond at Z around 0.025) and the VA84 so¬ 
lution (square with error bars) with the H80 models are 
shown. The Z value of VA84 (derived from the Nissen 
(1981) calibrations) is Z = 0.017 (vertical line). Two other 
points (filled triangle and filled diamond) are also shown 
at Z = 0.017, with the same ages as the (open triangle) 
and (open diamond) isochrones respectively. 

a remark has to be made. The masses of the TZ For sys¬ 
tem are known with great accuracy (better than 3% for 
the primary and almost 1% for the secondary) and provide 
stringent tests: we should derive the same masses with the 
theoretical tracks. Thus, if we consider TZ For B which is 
the strongest test and fix the mass of the Geneva tracks to 
its measured value. Fig. |T^ clearly shows that none of the 
tracks succeeds in fitting this star whatever the metallicity 
used (cf. left panel). Alternatively, if we fix the metallic¬ 
ity to its observed value (Andersen et al. 1991), the mass 
of the tracks have to be reduced by 5cr to reproduce the 
mass of TZ For B (see right panel). This shows that the 
Geneva models do predict a metallicity consistent with the 
spectroscopic value, but if the metallicity is hxed they do 
not predict the correct mass any more with the required 
precision. 

4.4.7. VV Pyx [19] 

The best isochrone solution found with the Padova tracks 
gives (logt, Z) = (8.721, 0.008) while the Geneva tracks 
predict logt = 8.753±0.04 also at the same metallic¬ 
ity. The Granada tracks give solutions in good agree¬ 
ment with Padova and Geneva: (logt, Z) = (8.76±0.04, 
O.OlOiO.OOl), see Fig. 0- It is unfortunate that there is 
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Fig. 16. TZ For system [18]. The hottest component of TZ For, TZ For B, is a star whose mass is known to 1% 
accuracy (1.949±0.022 Mq). Geneva tracks of several metallicities (with a mass equal to the TZ For B mass) are 
shown, but none of them fits this component {left panel). In order to fit TZ For B with the observed metallicity, the 
mass has to be reduced by 5cr {right panel). 


no metallicity indications, because Andersen, Clausen & 
Nordstrom (1984a) derived a solution with a metallicity 
twice as large (logt, Z) = (8.6, 0.015) using the H80 mod¬ 
els. The difference in the solutions certainly arises from 
the different physical ingredients of these tracks with re¬ 
spect to the sets we have used, as explained by Andersen 
et al. (1984). A measure of [Fe/H] would therefore be an 
ideal test here, even though spectra may be contaminated 
by the visual companion, a main—sequence A5—A7 star, 
itself perhaps a spectroscopic binary (Andersen, Clausen 
& Nordstrom 1984a). The results we derived from the 
HR diagram would be unchanged with the TeffS of Ribas 
et al. (2000), because the revision is tiny (ATeff=22 K). 
An attempt to derive [Fe/H] from the synthetic BaSeL 
photometry (Lastennet et al. 1999a) suggests a metal¬ 
licity markedly smaller than solar ([Fe/H]< —0.45, i.e. 
Z<0.007), but the fit of Stromgren colours and logg is 
bad. While this bad fit may reveal a problem of the BaSeL 
models, incorrect colours may be an alternative explana¬ 
tion, supported by the possible contamination of the visual 
companion. 


4.4.8. V1647 Sgr [22] 

Andersen & Gimenez (1985) used the H80 models to de¬ 
rive best fits between {\ogt, Z,Y) = (8., 0.03, 0.27) and 
(8.30, 0.02, 0.23). With the Geneva tracks, we get an age 
consistent with these results (logt ^ 8.24), but a solar 
metallicity {Z = 0.017), even though a larger metallicity 
(say Z = 0.023) at a slightly smaller age is not excluded 
by the y^-contours. The Padova models agree with these 
values (see Fig. logt < 8.40 and Z = 0.014to;oj2 
as well as the Granada tracks : log! £ [7.50, 8.62] and 
Z = 0.0161q qq 2 ' Andersen & Gimenez (1985) quote a pos¬ 
sible observed range for Z between 0.02 and 0.04, which 
is consistent with the 3 sets of tracks used here, although 
smaller values are permitted. If the lower limit of Z > 0.02 
is confirmed, the Granada tracks may not be able to ht 
the properties of this system. 

Note however that the rotation velocities of the com¬ 
ponents are not negligible (nsini)^ = 80 ± 5 km s“^ and 
{vsmi)B = 70 ± 5 km s“^ (Andersen & Gimenez 1985) 
and that isochrones incorporating rotational effects would 
give slightly younger ages and larger metallicities. It has 
also been pointed out that the parameters of this system 
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Fig. 15. PV Pup [12], Padova tracks. Isochrones for the 
best fit (filled circle at (logt, Z) = (5., 0.04)) and for points 
at Icr away (triangle), and for four coeval points (open cir¬ 
cles) with the best age solution but with different metal- 
licities {Z = 0.030, 0.035, 0.045 and 0.050). The VA84 
solution (square with error bars) with the H80 models is 
also shown, along with the Z = 0.017 (vertical line) de¬ 
rived from Stromgren indexes (VA84 with Nissen (1981) 
calibrations). Two other points (filled triangle and filled 
diamond) are also shown on the Z = 0.017 line. 

may have been perturbed by the presence of its bright vi¬ 
sual companion, and this has not been taken into account 
in the error budget (Ribas et al. 1998). 

4.4.9. YZ Cas [25] 

This system has been studied a lot previously, and Fig, 
and Table || summarize the results obtained so far. Except 
for the measures of Lacy (1981), all sets of models (what¬ 
ever the stellar parameters used) agree in supporting a 
solar metallicity (Z range between 0.015 and 0.020). As 
quoted in Tab. g a solar composition is supported by 
photometric determination on the B component, a F2V 
star (Ribas et ah, 2000). Even though its components are 
not very evolved, YZ Cas is a stringent test of the models 
because of the relatively large difference in its masses. 
The metal abundance determination of ROO is therefore 
a strong support to the quality of the tracks in the ~ 
1.3-2.3 Mq range. 

Apart from this very good agreement between models 
and data, we have to mention that there is still an open 
problem concerning the primary component, YZ Cas A, a 
metallic Am star. While a disagreement may be expected 



Fig. 17. VV Pyx system [19]. Comparison between the 
best fits obtained using the Geneva (circle), Padova 
(square) and Granada (triangle) tracks, superimposed on 
the confidence regions of the Geneva tracks. The best so¬ 
lutions of Andersen, Glausen & Nordstrom (1984a) us¬ 
ing the H80 models (diamond) and Pols et ah, 1997 (er¬ 
ror bars) are also shown. The photometric constraint on 
[Fe/H] (Lastennet et al. 1999a) is shown as a solid line 
with arrows. 

between the atmospheric metallicity (.^atmos.) and its in¬ 
trinsic composition (Z, as derived from the models) due 
to an enhancement of the metal abundance in its exter¬ 
nal layers, Lastennet, Valls-Gabaud & Jordi (2001) found 
some discrepancy between the Zatmos. determinations. 
While photometric methods suggest that Zatmos is solaiQ 
Zatmos derived from lUE observations by De Landtsheer & 
Mulder (1983) is still high (Zatmos.=0.036±0.005) despite 
the attempt of Lastennet, Valls-Gabaud & Jordi (2001) 
to revise some of the assumptions made by the lUE-based 
spectral study. More detailed analyses are needed to un¬ 
derstand this disagreement. 

4.4.10. (3 Aur [27] 

In this system we can compare the effects of increased 
precision in the measures of its parameters. Fig. ^ uses 
the logTeff and logL parameters listed by Andersen 
(1991) (upper panels) and the more recent one^ from 
Nordstrom & Johansen (1994b, hereafter NJ94b) in the 
lower panels. The former values give Z = 0.025 and logt= 
8.614 using the Geneva tracks, similar to the best fit val¬ 
ues found for each component separately^. These results 
are inconsistent with the constraints in metallicity given 
by Toy (1969) : —0.03 < [Fe/H] < 0, and the disagreement 

^atmos=0.017/Q qJy (Lastennet, Valls-Gabaud & Jordi, 
2001) and Zatmos<0.022 (L99a). 

Not modified by the TeS revision of ROO. 

(logt,Z)A = 8.587, 0.029 and (logt,Z)B = 8.640, 0.023. 
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Table 2. Comparison of age (logt, in yrs.), metallicity (Z), and helium abundance (T) determinations obtained for 
YZ Cas [25]. 


Tracks 

logt 

Z 

y(t) 

Reference 

Mengel et al. (1979) 

8.544T0.004 

0.027±0.003 

0.30±0.02 

(1) 

Hejlesen (1980a,b) 

Q ^1 q+U.U4U 

O.UlO_Q Q44 

0.03±0.01 

0.30±0.04 

(1) 

Mengel et al. (1979) 

8.8±0.1 

0.010T0.005 

0.30±0.05 

(2) 

Hejlesen (1980a,b) 

8.5T0.1 

0.015T0.005 

0.335±0.050 

(2) 

de Loore et al. (1977) 

8.65±0.10 

0.015±0.005 

0.335±0.050 

(2) 

Celikel & Eryurt-Ezer (1989) 

8.77±0.02 

0.021^ 

0.240^*^ 

(3) 

Pols et al. (1998) STD 

8.66±0.02 

0 m 4+0-004 
U.Ui+_0.003 

0.2681?:°°? 

(4) 

Pols et al. (1998) OVS 

8.70±0.02 

0 01 fi+0 006 

U.UiO_0.004 

n 979+0.012 

U.Z/^_0.008 

(4) 

Cg(“) 

8.654T0.014 

0.016±0.004 

0.270±0.026 

(5) 

CG92('’) 

Q 70+0.07 

o. ( 'J_o.i9 

o.oi7lH« 

U.Z,/ G)_Q QQg 

(6) 

Geneva 

8.651°;]? 

0.0181?:“? 

0.2931?;?]? 

(6) 

Padova 

8 . 62 in? 

U.UZU_o.o07 

0.2801?;??? 

(6) 

Stromgren photometry (B comp.) 


0.018-0.020 


(5) 


Claret (1995), Claret & Gimenez (1995), Claret (1997), Claret & Gimenez (1998); Claret & Gimenez (1992); Assumed composition. 
(1) Lacy (1981); (2) De Landtsheer & De Greve (1984); (3) Celikel & Eryurt-Ezer (1989); (4) Pols et al. (1997); (5) Ribas et al. (2000): 

[Fe/H]—0.03, i.e. 2^—0.018 (assuming 0.017, Grevesse 1997, priv. comm.) or Z—0.020 (assuming 2^0—0.0188, Schaller et al. 1992); (6) This 

work. Only the GG models allow an independent determination of Y and Z. Otherwise, Y is derived from a fixed law once Z is determined. 


is even worse with the Padova tracks, since they give best 
fit models a,t Z = 0.035 and logt = 8.537, although at the 
la level they may be compatible. 

The NJ94b analysis yields two components which are 
well separated in the HR diagram and the constraints in 
the (Z, logO diagram are much better, as the lower pan¬ 
els of Fig. Ily show. The uncertainty regions are greatly re¬ 
duced and the solutions are shifted towards smaller metal- 
licities. Nordstrom & Johansen (1994b) obtained Z = 
0.015 and an associated age of logf = 8. 750 The so¬ 
lutions obtained here with the Geneva {Z = O.Old^'^g QQ®), 
Padova {Z = O.OOOto'osgj see Fig. lower panels) and 
Granada {Z = 0.017lo oQg) tracks are all in good agree¬ 
ment and consistent with the metallicity constraint of Toy 
(1969). Note however that there is an island of solutions 
at high metallicity and smaller age that remains with the 
Padova tracks, giving rise to isochrones that are statis¬ 
tically indistinguishable from the lower Z ones (see the 
left bottom panel in Fig. ^4 )- This island is in principle 
excluded because the implied metallicity is much larger 
than the observed one, although a modern spectroscopic 
determination would be useful to confirm Toy’s (1969) val¬ 
ues. 

4.4.11. V1031 Ori [28] 

All the models converge towards the same solution for 
this system. Andersen, Clausen & Nordstrom (1990) used 
two sets of tracks and derived {logt, Z,Y) = (8.74, 0.024, 


As also noted by Nordstrom & Johansen (1994b), this sys¬ 
tem is about 8 times younger than the Sun, and yet has a 
similar metallicity. 


0.27) with VandenBerg’s (1985) models and (logt, Z, Y) = 
(8.85, 0.02, 0.28) with the older generation of Geneva mod¬ 
els (Maeder & Meynet 1988, 1989). The updated Geneva 
tracks give a best fit at (logt, ^) = (8.771, 0.021), while 
the Padova tracks indicate (logt, Z) = (8.747, 0.020), 
in agreement with Pols et al. (1997), (logt, Z) = (8.79, 
0.023), and Ribas et al. (2000), (logt, Z) = (8.846, 0.016). 
Granada models favour coeval solutions with a larger 
metallicity (Z ss 0.027), although at the la confidence 
level there is no disagreement with the other sets of tracks, 
and with the photometric constraint derived from the 
BaSeL models on V1031 Ori B (Z<0.029, Lastennet et 
al. 1999a). Nevertheless, another photometric determina¬ 
tion (Ribas et ah, 2000) suggests a sub-solar metallicity 
(Z~0.010) which would imply a systematic problem in all 
the above-mentioned theoretical models for main-sequence 
stars around 2.3-2.5 M©. The CG models used by Ribas et 
al. (2000) would be marginally consistent while the others 
are clearly ruled out. Once again, a detailed spectroscopic 
analysis would be needed before further inferences. 

4.5. Systems with at least one star more massive than 
3 Mq 

4.5.1. GG Lup [36] 

For this system we use the Andersen et al. (1993) data 
which update the A91 review. Even though the positions 
of the components in the HR diagram have not changed, 
the accuracy has increased two- or three-fold. This im¬ 
provement results in smaller uncertainty region^0 but 

The A91 parameters give Zceneva £ [0.011, 0.040], Zpadova 
G [0.010, 0.062] and Zoranada £ [0.010, 0.030], while the ones 









E. Lastennet & D. Valls-Gabaud: Detached binaries in the HR diagram 


21 




Fig. 18. V1647 Sgr system [22], Padova tracks. Upper 
panel: contours in the (metallicity, age) plane. Note the 

shift in metallicity between the best solution from Padova 
tracks (open circle at {logt,Z) = (8.38, 0.014)) and the 
best fit by Andersen & Gimenez (1985) (square with error 
bars) using the H80 models. The two other points (trian¬ 
gle and diamond) define the two corresponding isochrones 
(dashed line and dot-dashed line respectively) in the lower 
panel. Lower panel: Best fitting Padova isochrone (solid 
line), along with two isochrones at Icr. 


still highly elongated along the Z direction (Fig. ^l|) . 
Andersen et al. (1993) derived a best fit isochrone with 
(^, logt) = (0.015, 7.30) with the Granada modelsp^, 
claiming that this metallicity is in agreement with unpub¬ 
lished preliminary determinations by Clausen. Our results 
are consistent with these values at the Icr level, but a much 
wider range in Z is allowed by the evolutionary tracks. 
Basically all isochrones with ages younger than 5x10^ yr 
and a metallicity between Z=0.014 and Z=0.037 can fit 
reasonably well this system in the HR diagram. These re¬ 
sults are unchanged if we consider the recent Tea revision 
from ROO. On the other hand, Lastennet et al. (1999a) 
provide a slightly hotter Tes (by ~0.9%) and a smaller 
uncertainty (by a factor of 3) on the coolest component]^ 
which would imply smaller contours in Fig. excluding 
all the solutions with Z>0.025. 


from Andersen et al. (1993) produce Zoeneva £ [0.015, 0.037], 
ZPadova € [0.014, 0.035] and Zcranada € [0.015, 0.030]. 

P97 obtained similar solutions: Z=0.015±0.002, 
logt =7.18 )Iq 22- ROO obtained no solution in the (Z, Y) 
range covered by the CG models. 

log Tefr^d.OdblQ Q?!, to be compared with A91 or ROO; log 
Teft=4.041±0.024. 




Fig. 19. YZ Gas [25], Geneva tracks (open circle symbol 
for the best fit, and 1-, 2- 3-a contours). For comparison, 
the results of previous studies are indicated: Lacy (1981) 
using the MDSG79 tracks (filled square) and the H80 mod¬ 
els (filled triangle); De Landtsheer & De Greve (1984) 
with the MDSG79 (square), the H80 (triangle), and the 
De Loore et al. (1977) models (cross); Celikel & Eryurt- 
Ezer (1989) (filled circle). More recent studies from Pols et 
al. (1997) (open star), and Ribas et al. (2000) (diamond) 
are also shown. A metallicity estimate from photometry is 
indicated with vertical dotted lines. See Tablefor further 
details. 

4.5.2. V539 Ara [41] 

This system was studied by Andersen (1983), and more re¬ 
cently by Glausen (1996), allowing to study the influence 
of more accurate data on the best fit isochrones. Figure]^ 
uses the Andersen (1983) data in the top panels and the 
Glausen (1996) values in the lower panels. Both Geneva 
and Padova tracks give a smaller Z and a larger age than 
previous determinations^ although they all lie on the Icr 
contours. The updated values decrease somewhat the un¬ 
certainty area in the metallicity-age plane, but a large 
range of Z is still allowed (see Fig. ^), despite the two¬ 
fold increase in the accuracy of the data. Unfortunately, 
as indicated by Clausen (1996), attempts to measure the 
metallicity are hampered by the relatively large rotational 
velocities of its components, respectively 75 km s“^ (com¬ 
ponent A) and 48 km s“^ (component B). 

The best fit values obtained by Clausen (1996) using the 
Granada tracks are (logt, Z, Y) = (7.65, 0.02, 0.28), iden- 

e.g. Clausen (1979) with the Hejlesen et al. (1972) tracks, 
De Greve (1989) with or without overshooting (tracks from 
Prantzos et al. 1986). 
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Fig. 22. V539 Ara system [41]. Top panels: metallicity-age contours and HR diagram with the Andersen (1983) data. 
The best fitting isochrone with the current Geneva tracks (open circle, at (logt, Z) = (7.7, 0.01)) gives a smaller Z and 
larger age than previous determinations (square: Clausen 1979 [with the Hejlesen et al. 1972 tracks]; open triangle : De 
Greve 1989 with overshooting; filled triangle : De Greve 1989 without overshooting [the latter two with the Prantzos 
et al. 1986 tracks]). Lower panels: The same as above, but with the more recent values from Clausen (1996) and tracks 
from Geneva (G.) and Padova (P.). The diamond represents the best fit value obtained by Clausen (1996) using the 
Granada tracks. The open circles (G. at (logt,Z) = (7.58, 0.017) and P. at {logt,Z) = (7.62, 0.015)) give the best 
fits from Geneva (full line) and Padova (dashed), while typical Padova Icr isochrones (filled circles) are also shown for 
comparison (dot-dashed lines). 


tical to our results (see Table ||) , but there are many other 
possible solutions given approximately by the confidence 
regions {logt,Z) = (7.65to 2 i yr, 0.020to oJ8). This point 
emphasizes the power of the method used in the present 
work: the confidence regions are essential to properly as¬ 
sess and understand the results of the tests. 

4.6. Very massive stars: M > 10 Mq 

As quoted in 1987 by Hilditch & Bell, there were not many 
accurate systems in this mass range, and this is unfortu¬ 
nately still true: only 7 systems of our working sample fall 
in this category, namely (by increasing order of mass) GW 


Cep [43], AH Cep [51], V478 Cyg [44], Y Cyg [50], EM Car 
[45], V3903 Sgr [52] and DH Cep [53]^. The main reason 
is that the MS lifetime of such massive stars is short, and 
so these binaries are often observed during their interact¬ 
ing phase of evolution, excluding them from our sample of 
well-detached systems. 

All these systems have good fits, and the more massive 
ones (> 15 Mq, which then excludes CW Cep and AH 
Cep) see m to yield low metallicities. As mentioned earlier 
on (j |3.l| ), DH Cep may have overflown its Roche lobe, 
and we exclude it from further discussion (even though 

As quoted in the last one does not match the 1-2% 
level of accuracy of the core sample studied in this paper. 
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Fig. 20. [3 Aur [27], contours with Padova models. Upper 
panels : assuming the A91 values in the HR diagram, best 
fits with Padova (circle), Geneva (triangle) and Granada 
(open star) tracks. Note that the Granada value is at the 
upper limit available in this set {Z = 0.03). Lower panels 
: Using the updated parameters from NJ94b, the best 
fit, high metallicity, Padova isochrone (open circle) is con¬ 
sistent with the best fit isochrones from Geneva (trian¬ 
gle) and Granada models (open star), even though their 
metallicities differ by a factor of 3. Also shown is a Padova 
isochrone (dashed line) of solar metallicity (open square). 
The diamond represents the NJ94b solution (using the 
Granada tracks). The observed constraints on metallicity 
(Toy 1969) between [Fe/H] =—0.03 and [Fe/H] = 0 are 
indicated with vertical lines. 


it follows this trend of lower Z as well). Previous studies 
with different fitting methods also show the same trend: 
extrapolated solutions at Z = 0.009 for V478 Gyg, Y Cyg 
and EM Gar, with a possibly even lower Z for Y Cyg 
according to its bad fit (Pols et ah, 1997), Z =0.005 for 
EM Car and Z =0.010 for V3903 Sgr (Ribas et ah, 2000). 

Can new and carefully derived temperatures change 
this low metallicity trend? The revised temperatures of 
Ribas et al. (2000) leave unchanged the Teff of V478 Cyg, 
EM Car, Y Cyg and V3903 Sgr and decrease by only 0.003 
dex both components of GW Cep, therefore same results 
are expected with these new determinations. Obviously, 
results that can be derived for 0-B binaries in the HR di¬ 
agram have to be taken cautiously due to the more uncer¬ 
tain Teff-scale for massive stars, and so definitive conclu¬ 
sions are premature. Moreover, because metal abundance 
is difficult to determine in hot stars, it becomes difficult 
to discriminate between models. For this reason, EM Car 
deserves more discussion because it is the only system 


Fig. 21. GG Lup [36], Geneva tracks contours. Best fits 
with Geneva (circle), Padova (square) and Granada (tri¬ 
angle) tracks are shown, as well as the best fit by Andersen 
et al. (1993) (diamond, ACG93) with the Granada tracks. 
The data points in the HR diagram are from Andersen et 
al. (1993). 

with metallicity indications. Our theoretical predictions 
favour low metallicities: Z=0.003 (Geneva models) and 
Z=0.004 (Padova models), in agreement with Pols et al. 
1997 (Z=0.009, extrapolated solution), and Ribas et al. 
2000 (Z=0.005±0.002)Q 

The metallicity constraint on EM Car is [Fe/H] <0.10, 
as inferred from the photometry of its secondary com¬ 
ponent, and favours low metallicity solutions ([Fe/H] = 
—0.90, i.e. Z~0.004). This result is obtained from pho¬ 
tometric calibration (BaSeL library) matching the very 
accurate surface gravity (log g=3.928) and the observed 
dereddened colours (b—y)o, mp, and Cq (Lastennet et al. 
1999a), and supports the inferred low-Z trend mentioned 
above. 

Alternative explanations for low metallicities in mas¬ 
sive binaries may come from the fact that the models may 
be too dim, so low metallicties models can fit the data 
without a proper physical reason (Young et al. 2001). 
However, while Young et al. (2001) and our results do 
use the information on the luminosity, this do not explain 
the low-Z results of Pols et al. and Ribas et ah. Therefore 
the heavy element abundance may be actually small for 
the systems of the sample with masses larger than ~ 15 
Mq. More accurate data for massive stars are needed to 
solve this issue. 

However all these stellar evolution models might be sys¬ 
tematically wrong for very massive stars (rotation, mass loss 
rates, diffusion effects). 
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Fig. 23. Hya system [34], Geneva tracks. We use 3 different x^-estimators to fit both components : (1) a fit 
matching Teff and L (our standard formulation) (upper panels), (2) a simultaneous fit of Tea, L, M and R (middle 
panels), and (3) a fit involving only M and R (lower panels). For each x^ estimator the confidence regions in the 
metallicity-age plane are given (left panels), along with the resulting HR diagrams (middle panels) and Mass-Radius 
diagrams (right panels). For comparison, the minima obtained from each method (filled circle, square and triangle for 
methods (1), (2) and (3), respectively) are plotted on the upper left panel. 


5. General discussion 

We discuss here the implications of the careful determi¬ 
nations of age and metallicity that we have obtained for 
most of the systems of our sample. 


5.1. Possible systematic trends in age and metallicity 


Before attempting to draw general conclusions from these 
results, it is worth asking whether the inclusion of further 
information would improve the determinations of age and 
metallicity. After all, as we showed in Section § 3.3, when 
fitting a single isochrone to the system, the degeneracy 
area in the age-metallicity plane was in general reduced, 
due to the extra information added (two luminosities and 
two effective temperatures). Would the inclusion of mass 
and/or radius, into the x^ functional reduce the uncer¬ 


tainties in the final ages and metallicities? We compare 
here three different formulations for the x^ functional : 
(f) our standard expression, based only on Tea and L; (2) 
a functional minimizing the residuals using all the infor¬ 
mation available, namely the four quantities Tea, L, M 
and R] and (3) a functional using only the most direct 
quantities M and R. This last formulation allows us to 
check the predictive power of the theory on the other two 
parameters, Tea and L. As in ^3.3| , we will take the system 
X^ Hya [34] as an example. 

Each formulation is illustrated on Fig. with a Z-t 
diagram, a classical HR diagram and a M — R diagram. 
It is immediately clear that in the (Z, t) plane the con¬ 
fidence areas are roughly similar in shape and extension, 
although method (2) provides slightly tighter constraints, 
as expected, at the 3 (t level. However in this case there are 
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several “islands” at the Icr level, giving rise to multiple, 
unconnected solutions, which are difficult to understand 
given Vogt’s theorem. This multimodality also indicates 
that 1(7 solutions may be vastly underestimated if a lo¬ 
cal minimum (and not the absolute one) is found instead. 
These islands are even narrower when using M and R only 
(method 3). While method (3) makes a good fit to R and 
M, it fails to provide an excellent isochrone in the HR 
diagram. Likewise method (2), which has no predictive 
power in this sense, yields somewhat worse fits in the HR 
diagram than method (1) which is therefore our preferred 
one for providing the best possible fit in the HR diagram 
and for predicting properly masses and radii. 

At first sight, these results might be surprising: the ex¬ 
tra information added (method 2) does not yield a vast im¬ 
provement on the solution, nor do the more precise masses 
and radii (method 3). The addition of more information 
(method 2) imposes much tighter constraints on the mod¬ 
els, and simply reveals their deficiencies in a better way. 
There are two possible additional reasons for this. First, 
all these quantities are correlated, and hence bring less 
information than naively expected. When we added the 
two components in §3.3 fitting a single isochrone, the cor¬ 
relations between the parameters of one component and 
the parameters of the other remain, but are in general 
much smaller than those within a single component, and 
therefore yield better constraints. Another possible reason 
is that these systems remain close to the ZAMS, where 
there are well-defined mass-radius and mass-luminosity 
relations, implying even stronger correlations, and thus 
even less information. Systems with components in widely- 
separated stages of evolution would probably make much 
larger differences. With the current sample of binaries, 
however, these different formulations yield somewhat sim¬ 
ilar results, with method (1) being preferred for yielding 
the best fits to L and Tes, and predicting correctly M and 
R. 


An important motivation for this choice is that for single 
stars and less constrained systems such as visual binaries, 
one does not have M and R, whereas L and Tea can be 
determined in many cases. 

Another question that arises is whether the solutions 
found here are compatible with previous analyses, which 
used not only different sets of tracks, but also different 

formulations. By and large, the {Z,t) solutions are 
very similar and statistically consistent with each other. 
But are there any systematic trends between the different 
analyses? We compare here in detail our results with those 
obtained by Pols et al. (1997) and by Ribas et al. (2000). 
We do not take into account the results from Young et al. 
2001 because they fixed all the metallicities to the solar 
value. The sub-sample in common includes only 14 sys¬ 
tems mainly because we are restricted to Z<0.03, the up¬ 
per limit considered in P97 and ROO. We plot the relative 
differences (in the sense results with Padova, Granada, 
ROO or P97 minus results derived from the Geneva tracks) 
in age on Fig. Not only the agreement between the 5 
different estimations are very good, but there are no sys- 






Fig. 24. Relative errors in age (log t) derived with Padova 
models (a), Granada tracks (b), ROO (c) and P97 (d), as 
a function of the age predicted by the Geneva tracks. The 
dotted line is Alog t = 0 and the dot-dashed lines are linear 
regressions to the data points. 


tematic effects: the average relative error in age is —0.49%, 
0.32%, 0.29% and 0.35%, with a dispersion of 2.2%, 1.0%, 
1.3% and 0.5% for panels (a), (b), (c) and (d), respectively. 
This is a very good result for the theoretical models con¬ 
sidering the large individual error bars. 

The error in metallicity (Fig. 25) is much larger, but 
results are consistent within 20% for most cases, even 
though uncertainties remain large. The metallicity pro¬ 
vided by the Geneva tracks is compatible with the one 
given by the Padova tracks (Fig. p^a). There is a slight 
systematic trend of the Geneva (and hence Padova) tracks 
with the other ones (Granada, P97 and ROO) in the sense 
that the low metallicity systems are given a smaller metal¬ 
licity by the Geneva (and Padova tracks). However this 
trend is produced by two stars only {( Phe [35] and EM 
Car [45]) and seems hardly significant because the P97 
and ROO solutions are extrapolated for such low Z-values. 


5.2. Metallicity constraints 

Given that all these independent methods seem to predict 
consistent ages, how do the metallicity values compare to 
the observed ones? Obviously the most discriminant sys¬ 
tems are those, as pointed out by Andersen (1991) and ear¬ 
lier references therein, which have relative errors smaller 
than 2% in mass, less than 1% in radius, less than 2% 
for temperature and less than 25% for [Fe/H]. Nearly all 
the binaries presented in this paper verify the two first 
criteria, but only 32 individual stars reach the target ac¬ 
curacy on Teff. In addition, we want both components in 









































26 


E. Lastennet & D. Valls-Gabaud: Detached binaries in the HR diagram 






Fig. 25. Error in metallicity derived with Padova models 
(a), Granada tracks (b), ROO (c) and P97 (d), as a function 
of the metallicity predicted by the Geneva tracks. The 
dotted line is AZ = 0 and the dot-dashed lines are linear 
regressions to the data points. 


a system to be within these limit of a 2% error in Tgff, so 
we have to exclude 4 systems ([18], [35], [46], [59]). The 
revision proposed by Ribas et al. (2000) affects some of 
these stars, so for extra precaution we also exclude the 
systems revised by more than 2% ([7], [20], [24] and [33]). 
Among this reduced subsample of 10 systems^, we have 
to further exclude systems showing a bad fit in the HR 
diagram (see p.3| ) or those for which models are unable 
to predict correctly the mass and/or radius (^.4|). The 7 
systems selected ([3], [6], [8], [9], [11], [50], [56]) provide 
the best possible tests for the comparison of [Fe/H]. 

Many attempts have been done in the past to measure 
[Fe/H] in double-lined eclipsing binaries. A summary is 
presented in Table |^. It is somewhat distressing to note 
that from all the systems listed in Table only UX Men 
[6] and DM Vir [8] belong to the sample selected above 
with strict accuracy criteria. 

An inspection of Table ^ shows that the derived and 
observed metallicities are not always in perfect agreement 
but are consistent within the quoted error bars. The most 
significant discrepancy appears for the system UX Men [^, 
where the theoretical models predict too large Z-value^, 
even if they are consistent at the 1-cr level as shown in 

[2], [3], [4], [6], [8], [9], [11], [50], [56] and [57]. 

Even if Ribas et al. (2000) mention that their solution is un¬ 
reliable (probably due to the Z=0.03 upper limit of the models 
they used), it is worth noting that they also obtain an extrap¬ 
olated large metallicity (Z=0.039). This is the result we would 
have obtained assuming the lower Tetrs they adopt. 


Fig. 0. Whilst the disagreement may well come from the 
physical ingredients assumed in the tracks, we note that 
larger effective temperatures (by ~100-200 K) for each 
component would give a good agreement. Assuming the 
[Fe/H] value from spectroscopy (Andersen et al., 1989) 
and E(b—y)=0.07 (the colour excess derived by Lastennet, 
Cuisinier & Lejeune, 2002), the (b—y)-Teff calibration 
from the BaSeL models points to an increase by ~250 
K and 200 K for UX Men A and B respectively^. It is 
then likely that the most significant discrepancy shown in 
Table might be accounted for by the Tgff scaling. 

Due to the intrinsic difficulty to derive metallicities 
from spectroscopic analyses, photometric calibrations can 
be very helpful. Lastennet et al. (1999a) presented some 
metallicity-dependent temperatures determinations of 20 
EBs from Stromgren synthetic photometry (BaSeL mod¬ 
els), and Table 5 of Ribas et al. (2000) give determinations 
from other photometric and spectrophotometric methods. 
We exclude from these studies the results based on compo¬ 
nents with atmospheric anomalies (mainly Am stars) and 
present a summary in Table There are only 3 systems 
in common: WX Gep [30], TZ Men [29] and V1031 Ori 
[28]. The agreement is good for [28] and [29] (mainly be¬ 
cause the BaSeL results are weakly constraining for these 
systems), but not for WX Cep [30]. However, inspection 
of Fig. 4 of Lastennet et al. (1999a) shows that the BaSeL 
models are consistent with the ROO result at the 3 -ct level 
for this system. 

Another special case is the binary GZ CMa [23] whose 
primary component is an Am star. As we mentioned in 
§3.2, the difference between the revised ROO TgffS and 
the A91 TgffS which we adopt for this system is larger 
than 4%, so we consider our results ([23] in Tab. ||) to 
be less reliable. This gives rise to different interpretations 
even if our results are marginally consistent with their 
lower limit: we predict (Tab. ^ basically solar abundance 
(Z~0.018), while Ribas et al. obtain a metal-rich solution 
(Z~0.032±0.006). As shown in Table both the ROO and 
L99a TgffS are cooler than the values of A91. Adoptiim 
these TgffS would give a more metal-rich solution in Tab. ^ 


Table 5. Comparison of Tgff determinations for GZ CMa 
[23]. 




log Teft (K) 


GZ GMa 

A91 

L99a 

ROO 

A 

3.945±0.017 

3.928±0.008 

3.927±0.017 

B 

3.931±0.017 

3.922T0.005 

3.914F0.017 


The photometric calibration used by Ribas et al. 
strongly supports a large value for the metallicity. Again, 

These results are sensitive to the reddening. Assnming 
E(b—y)=0.04 would give an increase of about 70 K for both 
components. 
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the temperature scale proves to be a critical issue, and yet 
this comparison may not be entirely relevant because the 
photometrically-derived Z is inferred from the primary 
component which is an Am star. An estimation from the 
secondary component would be far more secure, and the 
BaSeL models (Tab. indicate a strict upper limit of 
Z <0.017. There is clearly an anomaly: on the one hand, 
a metal-rich solution is obtained when adopting the most 
recent revised Tefjs, and on the other hand photometric es¬ 
timates from the secondary component yield a solar value. 

This discussion is relevant for the problem of the chem¬ 
ical enrichment ratio AY/AZ carefully derived by Ribas 
et al. (2000), because GZ CMa is the only binary of their 
sample at ^>0.026 (cf. their Figure 7). Since we suggest 
to shift the metallicity of GZ GMa to solar values (or even 
slightly less), the data points shown in their Fig. 7 would 
cluster even more. Although this would not change much 
their linear relationship between initial metal and helium 
abundances, the Z range would be much reduced and the 
results would be less constraining for chemical evolution. 
A detailed spectroscopic analysis of GZ GMa, and par¬ 
ticularly of the secondary component, a (normal) A-type 
star is much needed to resolve this issue. 


5.3. On the relative and absolute dating of stars 


Since the ages derived by the analyses done here seem 
to agree with each other and with other determinations 
(j jS.ll ) and that, in addition, the predicted metalli citie s 
are roughly consistent with their observed values (§ |5.2|) , 
we may consider that at least the relative dating of these 
stars is well rooted. 


This has reassuring consequences for the chemical evo¬ 
lution of the Galactic disc. The detailed study made by 
Edvardsson et al. (1993) presented precise photometric 
and spectroscopic data for a sample of 189 field stars, and 
showed that there is no tight age-metallicity relation. This 
was recently confirmed with a larger sample (Feltzing et 
al. 2001): at a given age, the scatter in metallicity is very 
large. 


In this paper we derived simultaneous metallicity and 
age estimates (listed in Table ||) for 60 eclipsing binaries, 
most of them being representative of single stars. Of course 
this sample is highly biased and cannot be used for chem¬ 
ical evolution studies, but it is interesting to see, given 
the very precise dating that is achieved by this analy¬ 
sis, whether these stars follow the same trend. The ages 
derived by Edvardsson et al. from isochrone fitting range 
from about r = log t = 9.2 to 10.25. We have no stars older 
than logt =10.25, the oldest being ~9.5 Gyrs old (HS Aur 
[2]), but many more younger than their limit. After ex¬ 
cluding stars with weak constraints on their age and/or Z, 
well-fitted stars with ages between 10® and 10® yrs present 
a spread in Z between 0.05 and 0.01, i.e [Fe/H] between 
~-|-0.49 and —0.23 (or even below if the low metallicity 
of VV Pyx [19] is real, see § 1.4.7 ). This trend suggests 


that the dispersion observed by Edvardsson et al. may 
well continue at younger ages. 

6. Conclusions 

We have developed a method to derive simultaneously 
ages and metallicities from binary stars and single stars. 
Instead of giving one solution (i.e. one isochrone defined 
by a unique Z-t pair), this method provides a set of pos¬ 
sible solutions, as defined by y® contours. Within them, 
Vogt’s theorem is (operationally) violated: stars of dif¬ 
ferent masses and chemical compositions lie on the same 
position in the HR diagram. We have applied this method 
to the largest sample of well-detached binary systems for 
which most stellar parameters are extremely well known, 
to test the power of the method which could be applied to 
samples with looser constraints. In addition the method 
allows to test the predictive power of stellar evolution the¬ 
ory, by taking only two of the well-constrained quantities, 
and predicting the others. The choice of luminosity and ef¬ 
fective temperature as primary quantities turned out to be 
the best one, with accurate determination of the param¬ 
eters of the isochrone and also with excellent predictions 
for masses and radii (dispersion in their relative errors 
smaller than 2% and 4%, respectively). The validation of 
the method through these systems makes it possible to 
apply it to other, less constrained systems such as visual 
binaries, for which catalogues are growing much faster. 

The critical issue of systematic errors due to the use 
of a particular set of evolutionary tracks turned out not 
to be very important, at least for the sample used here. 
While the absolute results can be different in individual 
cases, all three sets of tracks used here (Geneva, Padova 
and Granada) yield - in a statistical sense - consistent 
values and similar age-metallicity degeneracies. 

We have made a detailed comparison of our results 
with previous studies whenever possible. The main 
inconsistencies found and discussed in this paper are the 
following: 

a) We have shown that the TeffS of YY Gem [la] may still 
need to be revised to slightly cooler temperatures. The 
Brocato et al. (1998) models are not able to fit the GM 
Dra [lb] components (both masses below 0.25 Mq) in the 
HR diagram for the expected solar metallicity, while the 
Baraffe et al. (1995) models succeed. 

b) We obtain (bad) fits which suggest large metallicities 
for 4 systems which have at least one component in 
the 0.7-1.1 Mq mass range. This supports the results 
found by Popper (1997) and Glausen et al. (1999) 
which indicated this peculiarity. Whilst some tentative 
explanations are discussed (Teff-scaling, activity, mass 
transfer, etc...), it appears that no Geneva isochrone 
can fit simultaneously both components of HS Aur [2] 
in the HR diagram. The best fit obtained with Padova 
tracks yields a (questionable) metal rich solution, about 
twice solar. This provides an ideal observational test to 
disentangle these sets of tracks. 

c) If the metal content of DM Vir [8] is confirmed to 
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be similar to the one of the Hyades, the Padova models 
would fail to fit the system, while the Geneva and 
Granada ones would succeed. This system is another test 
case to check stellar evolution theory in stars of masses 
of about 1.4 Mq. 

d) For RZ Cha [10], we have updated the not very 
constraining [Fe/H] value obtained by Jprgensen & 
Gyldenkerne (1975) (the Hyades metallicity they adopted 
was an overestimation), giving a very good agreement 
with the Geneva tracks, and a good agreement with the 
Granada and Padova tracks. 

e) The three sets of theoretical tracks suggest a metal rich 
solution for the binary PV Pup [12] (~ 1.55 Mq for both 
components) which is in disagreement with indications 
of a solar metallicity (Z~0.017, Lastennet et al. 1999a), 
even if the solar metallicity remains possible at the 2-a 
level. The new Tes determinations of Lastennet et al. 
(1999a) provides a slightly better agreement with solar 
metallicity values. 

f) For the ~2 Mq evolved components of TZ For [18], the 
agreement of a metallicity around the spectroscopic value 
(nearly solar) is satisfactory for the Geneva tracks, but 
seems difficult to reconcile with the Padova models. In 
spite of the better agreement obtained with the Geneva 
models, we have shown that if the mass of the tracks is 
hxed to the measured value of TZ For B, none of the 
tracks succeeds in fitting this star, independently of the 
metallicity used. Alternatively, if the metallicity is hxed 
to its spectroscopic value, the mass of the tracks must be 
reduced by 5a (0.11 M©) to reproduce the mass of TZ 
For B. 

g) If the lower limit of a solar metallicity is conhrmed 
observationally for V1647 Sgr [22], the Granada tracks 
may not be able to ht the system, although the effects of 
rotation must be assessed hrst for an unambiguous test. 

h) The disagreement between determinations of metallic¬ 
ity in GZ CMa [23] may bear some consequences on the 
chemical enrichment ratio AY/AZ carefully derived by 
Ribas et al. (2000). While they hnd that GZ GMa is the 
more metal rich binary of their sample (Z=0.032), the 
BaSeL models suggest a solar or sub-solar composition 
(L99a). 

i) Our analysis predict that the heavy element abundance 
may be small (Z<0.010) for systems with components 
of masses larger than ~ 15 M©. The photometrically 
derived value for EM Gar [45] supports this pattern. 

We summarize the general results obtained as follows. 

1. There is a wide age-metallicity degeneracy in the 
htting of double-lined eclipsing binaries, as well as for 
their individual components and therefore also for stars 
in the held, no matter how accurate the measures are. 

2. This degeneracy limits, in practice, the empirical 
validity of Vogt’s theorem, assuming its premises give 
sufficient and necessary conditions for the existence of 
unique solutions to the equations of stellar structure. 

3. The systematic errors introduced by the use of a 
given set of tracks have been assessed by comparing the 


results obtained for three different sets (Geneva, Padova, 
Granada). The systematics are rather small, and in 
most cases the solutions are statistically compatible at 
the Icr level. This may be due, in part, to their similar 
metallicity-correlated He abundance. 

4. In spite of the high precision measures that lead to 
very small errors in the astrophysical parameters of these 
systems, the lack of metallicity determinations limits a 
precise measure (O.I dex) of the ages, except in a few 
odd cases. Since photometric determinations provide 
useful constraints for only few systems, the spectroscopic 
measure of the metallicity in the components of these 
systems is essential for this technique to be an efficient 
test of stellar evolution. 

Another way to bypass the difficulty to determine 
the metallicity of a particular binary would be to select 
binaries which belong to star clusters with a very well 
measured metallicity (as proposed by Glausen & Gimenez 
1996; see Lastennet et al. 1999b for an application to 
Hyades binaries, and Lastennet, Valls-Gabaud, Oblak 
2000). This is a promising area given the difficulties to 
derive accurate metallicities from detailed spectroscopic 
analyses of some eclipsing binaries. 
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Table 3. Summary of the observational metallicities from spectroscopic analysis (sp.) or photometric calibrations in 
the Stromgren (S) or Geneva (G) systems. As in Table T (col. 7) refers to the tracks used to derive the values Ztracks 
and [Fe/H]tracfes given in the 2 last columns. The high metallicity results {Ztracks) for V1143 Gyg [7] are probably 
unreliable: applying the revised Tea from ROO would give a smaller value {Z in the range 0.02-0.03). 
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U.Ui / _o .004 

+0.00l°;f2 

V1143 Gyg [7] 

-b0.08 


S 

A87 

0.020 

1 

o.o4oi°:°°^ 

+0.38^:“ 






3 

0.0501°:°°^ 

+0.491°;“^ 







4 

u.uou_o.oio 

+0.25^:;° 

TZ For [18] 

-b0.1±0.2 


sp. 

A91b 

0 021 

U.UOi_Q 007 

1 

0.015lHoi 

-0.05l°i° 


+0.14±0.15 


sp. 

A91b 

u.uoo_ 0.006 

3 

0.0381°°?^ 

+0.37l°“/y 


+0.1±0.15 

0 021 ^ 
U.U^l_Q QQ0 


A91b 

u.uz,T_o 006 

4 

U.UiO_o.o02 

-O-OStVol 

UX Men [6] 

-0.15 


S 

GG76 

0.012 

1 

0.031±0.009 

+0.27l°;l^ 


-b0.04±0.10 

0.019±0.004 9 

sp. 

A89 

0.018+°;°Si 

2 

0.032±0.008 

+0.29+°;l° 


-0.05±0.15 


S 

A89 

U.U10_o.o04 

3 

U.U04_o_o09 

+0.311°;!^ 


-0.04 


G 

KN90 

0.015 

4 

u.uz / _o.o07 

-rU.zi_o 04 


0.1 

0.024 


T97 





AI Phe [5] 

-b0.17±0.20 


S 

VH85 

u.uzo_o.oo9 

1 

0.010tHo2 

fl 9Q+0-04 

U.ZO_o.io 


-0.14±0.1 

0.012±0.003 

sp. 

A88 

U.U1Z_0.002 

2 

0.009t°°°t 

-0.28l°;J® 


-0.1 

0.012 


T97 

0.014 

3 

0.008tHoo 

-0.331°:),^ 







4 

O.OlliO.OOl 

-0.19±0.04 

PV Pup [12] 

-0.06 

0.017 

S 

VA84 

0.015 

1 

> 0.024 

> 0.16 






3 

0.040t°:°f2 

+0.391°;?® 







4 

> 0.021 

> 0.10 

V1647 Sgr [22] 


0.02-0.04 

7 

AG85 

0.07-0.39 

1 

U.Ui / _0.003 

+0.00l°;J® 






3 

-^^-0.002 

U.Uy_Q QQ 







4 

0.016tHo2 

—0 03+°-^*^ 
U.UO_o.o6 

DM Vir [8] 

0 .12±0.12 

0.023±0.006 

s 

A84b 

u.uz,z,_o 005 

1 

0.034t°;°?i 

+0.311°?? 


0.13±0.05 


s 

this work^ 

u.uz,o_o.oo3 

3 

0.04lin?3 

+o.4ol°:?° 






4 

0.030t°:°?? 

+0.251°;°° 


“ Values of Z as given in the papers of column 5 (Ref.). ^ Adopting Zq = 0.017 and using Eq. ^ 1 ct confidence interval, 

fitting both components except for TZ For [18]: fitting only TZ For B with the CG92 tracks (T=4). [Fe/H] = 0.00 for the 
binary and [Fe/H] = — 0.03 for the secondary component. Me Donald observations (A91b). ESO observations (A91b). 
Adopted mean by A91b. ^ Assuming Zq = 0.0169 (VandenBerg 1985). ® Assuming Zq = 0.0189 (Maeder & Meynet 1988, 
1989). ^ Updade inferred from [Fe/H] = [Fe/H]Hyades — [Fe/H]Ami, by adopting [Fe/H]Hyades=0.14±0.05 (Perryman et ah, 
1998). 
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Table 4. New Z determinations from photometric calibrations: BaSeL models (Lastennet et al. 1999a) and ROO 
(Ribas et al. 2000), both based on the Stromgren photometry of individual stars. The Z-values are derived according 
to the correction presented in §2.2. We exclude from these studies the results based on components with atmospheric 
anomalies (mainly Am stars). The column “Notes” gives the component(s) used (A for primary, B for secondary 
component). 


System 



BaSeL 



ROO 




[Fe/H] 


Note 

[Fe/H] 

z t 

Note 

BW Aqr 

[9] 

-0.20T0.25 

0 01 1+0.008 
U.U11_Q QQ4 

A & B 




AR Aur 

[48] 

> +0.00 

>0.017 

A 




GZ CMa 

[23] 

< +0.00 

<0.017 

B 




EM Car 

[45] 

< +0.10 

<0.021 

B 




YZ Gas 

[25] 

no constraints 


B 

+0.03 

0.018 

B 

WX Cep“ 

[30] 

> +0.05 

>0.019 

A 

-0.03 

0.016 

A & B 

GW Cep 

[43] 

< -0.15 

<0.012 

A 




RZ Cha 

[10] 

[-0.50,0.00] 

[0.006,0.017] 

A & B 




V380 Cyg*” 





-0.03 

0.016 


HS Hya'’ 





-0.17 

0.012 


KW Hya 

[15] 

+0.00+0.10 

0.017+0.004 

B 




GG Lap 

[36] 

< -0.20 

<0.011 

A 




TZ Men 

[29] 

-0.20+0.40 

0 01 1 +0-016 
U.Uii_o.oo6 

A 

-0.03 

0.016 

B 

V451 Oph 

[31] 

< +0.10 

<0.021 

B 




V1031 Ori 

[28] 

< +0.25 

<0.029 

B 

-0.25 

0.010 

A & B 

IQ Per 

[33] 

< +0.25 

<0.029 

B 




Al Phe 

[5] 

[-0.35,-0.15] 

[0.008,0.012] 

B 




C Phe 

[35] 

< +0.20 

<0.027 

B 




PV Pup 

[12] 

< -0.40 

<0.008 

A 





^ This work, assuming Ze=0 .017 and X© =0.713 (Grevesse, 1997, priv. comm.). “ According to L99a (cf. their Tab. 2 and Fig. 
3). Spectrophotometric determination, see Ribas et al. (2000) and references therein. 
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Table 6. Ages and metallicities for the best fit isochrones to the sample of 58 double-lined eclipsing binaries for the 
different sets of tracks used^. Systems YY Gem [la] and CM Dra [lb] are discussed in §4.2, T refers to the tracks 
used to derive the fit : (1) Geneva with overshooting, (2) Geneva without overshooting, (3) Padova, and (4) CG92. 
Maximum and minimum values are given at 68% c.l. (±1 a) for the fits of individual components, while the values for 
the combined system are at ±na. Ages consistent with the ZAMS are indicated by ’ZA’, that is, t = tzAMS in Eq- 0- 


# 

Name 


Component A 

Component B 

System A-l-B 

T 

T 

^min 

^maix 

Z 

^min 

■^max 

T 

^min 

^max 

Z 

^min 

■^max 

T 

^min 

^max 


^min 

^max. 

n 

2 

HS Aur 

1 

10.06 

9.95 

10.15 

0.040 

0.020 

0.040 

9.86 

9.24 

10.10 

0.040 

0.017 

0.040 

9.98 

9.94 

10.02 

0.040 

0.035 

0.040 

3 



2 

10.13 

9.93 

10.19 

0.040 

0.020 

0.040 

9.96 

9.28 

10.06 

0.039 

0.017 

0.040 

9.98 

9.92 

10.02 

0.040 

0.035 

0.040 

3 



3 

10.05 

7.43 

10.14 

0.049 

0.022 

0.100 

8.64 

ZA 

10.05 

0.073 

0.016 

0.100 

9.98 

9.87 

10.06 

0.049 

0.021 

0.050 

3 

3 

EW Ori 

1 

8.39 

ZA 

9.15 

0.040 

0.024 

0.040 

8.49 

ZA 

9.22 

0.040 

0.028 

0.040 

8.42 

ZA 

8.79 

0.040 

0.036 

0.040 

1 



2 

7.35 

ZA 

9.11 

0.040 

0.024 

0.040 

7.27 

ZA 

9.15 

0.040 

0.025 

0.040 

7.30 

ZA 

8.75 

0.040 

0.032 

0.040 

1 



3 

ZA 

ZA 

9.00 

0.051 

0.023 

0.069 

ZA 

ZA 

9.14 

0.053 

0.024 

0.079 

ZA 

ZA 

8.80 

0.051 

0.032 

0.065 

1 



4 

5.91 

ZA 

9.10 

0.030 

0.017 

0.030 

5.91 

ZA 

9.16 

0.030 

0.018 

0.030 

9.60 

ZA 

9.60 

0.030 

0.021 

0.030 

1 

4 

FL Lyr 

1 

9.36 

8.77 

9.55 

0.040 

0.020 

0.040 

9.43 

ZA 

9.75 

0.040 

0.021 

0.040 

9.38 

9.19 

9.50 

0.040 

0.027 

0.040 

1 



2 

9.37 

8.67 

9.55 

0.040 

0.021 

0.040 

9.55 

ZA 

9.78 

0.040 

0.018 

0.040 

9.40 

9.09 

9.50 

0.040 

0.028 

0.040 

1 



3 

7.73 

ZA 

9.50 

0.062 

0.020 

0.086 

3.79 

ZA 

9.93 

0.081 

0.029 

0.100 

7.69 

ZA 

9.41 

0.065 

0.021 

0.084 

1 

5 

AI Phe 

1 

9.66 

9.64 

9.76 

0.009 

0.009 

0.028 

9.60 

9.55 

9.68 

0.008 

0.004 

0.012 

9.66 

9.64 

9.67 

0.010 

0.008 

0.011 

1 



2 

9.71 

9.60 

9.73 

0.025 

0.008 

0.040 

9.60 

9.59 

9.66 

0.006 

0.006 

0.013 

9.64 

9.61 

9.65 

0.009 

0.008 

0.013 

1 



3 

9.68 

9.67 

9.70 

0.027 

0.022 

0.035 

9.60 

9.54 

9.70 

0.007 

0.006 

0.013 

9.60 

9.59 

9.63 

0.008 

0.008 

0.011 

1 



4 

9.67 

9.64 

9.68 

0.011 

0.010 

0.015 

9.65 

9.63 

9.68 

0.010 

0.010 

0.013 

9.66 

9.64 

9.67 

0.011 

0.010 

0.012 

1 

6 

UX Men 

1 

9.42 

8.90 

9.49 

0.033 

0.021 

0.060 

9.44 

8.80 

9.49 

0.030 

0.018 

0.062 

9.45 

8.95 

9.56 

0.031 

0.022 

0.056 

1 



2 

9.43 

9.15 

9.55 

0.031 

0.021 

0.040 

9.47 

9.10 

9.58 

0.027 

0.019 

0.040 

9.45 

9.29 

9.52 

0.032 

0.024 

0.040 

1 



3 

8.83 

6.86 

9.20 

0.060 

0.024 

0.092 

7.97 

7.33 

9.33 

0.069 

0.020 

0.087 

9.26 

9.15 

9.38 

0.034 

0.025 

0.044 

1 



4 

9.38 

9.08 

9.57 

0.030 

0.017 

0.030 

9.44 

9.12 

9.60 

0.024 

0.016 

0.030 

9.40 

9.27 

9.51 

0.027 

0.020 

0.030 

1 


^ The full table is only available in electronic form. The sample printed here is for orientation on contents and format only. 



















